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Abstract 

We review recent progress in the study of the anomalous baryon number non- conservation 
at high temperatures and in high energy collisions. Recent results on high temperature phase 
transitions are described, and applications to electroweak baryogenesis are considered. The 
current status of the problem of electroweak instanton-like processes at high energies is outlined. 
This paper is written on the occasion of Sakharov's 75th anniversary and will appear in the 
memorial volume of Uspekhi (Usp. Fiz. Nauk, volume 166, No 5, May 1996). 



CERN-TH/96-13 
INR-913/96 
March 1, 1996 



1 rubakov@ms2. inr.ac.ru 
2 mshaposh@nxth04.cern.ch 



Contents 

I Introduction 



2 Basics of anomalous non-conservation of fermion quantum numbers 



3 Baryon asymmetry: preliminaries 



4 Sphaleron rate at finite temperatures 

4.1 Qualitative discussion 

4.2 The Green's function approach 

4.3 The relation to the "probability flux" formulae 

4.4 Quantum versus classical rate 

4.5 The sphaleron rate in the broken phase .... 

4.6 Real time numerical simulations 

4.7 Strong sphalerons 

4.8 Concluding remarks 



5 Phase transitions in gauge theories 

5.1 Equilibrium approximation 

5.2 Simple estimates] 

5.3 The infrared problem and factorization 

5.4 Phase structure of 3d theory 

5.5 The lattice formulation 

5.6 Some lattice results 

5.7 Dynamics of the phase transition . . . 



6 Survival of primordial baryon asymmetry 



7 Electroweak baryogenesis] 

7. 1 Strength of the phase transition 

7.2 Sources of CP- violation in the EW theory and its extensions 

7.3 EW baryogenesis: how to state the problem 

7 A Uniform scalar fields 

7.5 Asymmetry from fermion-domain wall interactions 

7.6 Strength of CP violation and baryon asymmetry 



8 Instanton-hke processes in high energy collisions. 

571 Summary of perturbative analysis about the instanton 

3.2 Unitarity| 

3.3 From many — » many to few — > many 

S Conclusion 



1 



1 Introduction 



In his famous paper ffl] Sakharov discussed for the first time the possibility of explaining the 
charge asymmetry of the Universe in terms of particle theory. The paper was submitted to JETP 
Letters in September 1966, two years after the discovery of CP- violation in K° decays |2| and 
one year after the microwave black-body radiation, predicted by the Big Bang theory ||, was 
found experimentally |||. To explain baryon asymmetry, Sakharov proposed an approximate 
character for the baryon conservation law, i.e., baryon number non-conservation and proton 
decay. Three years later Kuzmin published a paper || where a different model leading to the 
baryon asymmetry was constructed. One of its consequences was another process with B non- 
conservation, namely, neutron-antineutron oscillations. Since that time the idea that baryon 
number may not be exactly conserved in Nature has been elaborated upon considerably, both 
in the context of the generation of the baryon asymmetry of the Universe 0, [7[ |], [| [10], |ll 



(for reviews see refs. [|12|, [13], [L4[]) and because of theoretical developments that have lead to a 
unified picture of fundamental interactions. In the mid-70's, grand unified theories with inherent 
violation of baryon number were put forward [15|, |16], [17|, 18, 19]. Almost at the same time it 



was realized [EG, pTT] that non-perturbative effects related to instantons [22] and the complex 



structure of gauge theory vacuum pOJ [23|, |24|] lead to the non-conservation of baryon number 
even in the electroweak theory; it has been understood later |25j that similar effects are relevant 
for the baryon asymmetry. 

In his paper []I| Sakharov writes: "According to our hypothesis, the occurrence of C asymme- 
try is the consequence of violation of CP invariance in the nonstationary expansion of the hot 
Universe during the superdense stage, as manifest in the difference between the partial proba- 
bilities of the charge-conjugate reactions." Today, this short extract is usually dubbed as three 
necessary Sakharov's conditions for baryon asymmetry generation from the initial charge sym- 
metric state in the hot Universe, namely: 

(i) Baryon number non- conservation. 

(ii) C and CP violation. 

(iii) Deviations from thermal equilibrium. 
All three conditions are easily understood. 

(i) If baryon number were conserved, and the initial baryonic charge of the Universe were 
zero, the Universe today would be symmetric, rather than asymmetric The statement of the 
necessity of the baryon number non- conservation was quite revolutionary at that time. Today 
it is very natural theoretically; still, lacking positive results from experiments searching for B 
non-conservation, the baryon asymmetry of the Universe is a unique observational evidence in 
favour of it. 

(ii) If C or CP were conserved, then the rate of reactions with particles would be the same 
as the rate of reactions with antiparticles. If the initial state of the Universe was C- or CP- 



3 Of course, there is a loop-hole in this argument, which Sakharov knew. The Universe may be globally 
symmetric, but locally asymmetric, with the size of the baryonic cluster of matter large enough (say, of the order 
of the present horizon size). The inflationary models of the Universe expansion, together with specific models of 
particle interactions may provide a mechanism of the local asymmetry generation, keeping the conservation of 
the baryon number intact fl3fl. 
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symmetric, then no charge asymmetry could develop from it 0. In more formal language, this 
follows from the fact that if the initial density matrix of the system po commutes with C- or 
CP-operations, and the Hamiltonian of the system is C- or CP-invariant, then at any time the 
density matrix p(t) is C- or CP-invariant, so that the average of any C- or CP-odd operator is 
zero. 

(iii) Thermal equilibrium means that the system is stationary (no time dependence at all). 
Hence, if the initial baryon number is zero, then it is zero forever. 

Clearly, the issue of the baryon asymmetry generation requires the development of many 
different areas of theoretical physics, such as model building, study of perturbative and non- 
perturbative effects leading to B-violation, finite temperature field theory and non-equilibrium 
statistical mechanics, theory of phase transitions. 

In this paper we do not aim to give a complete review of various theories of baryogenesis 
proposed so far. The reader may consult with a number of reviews on this subject fli2| , [T3| , 



6| ^7|> EI> HI' EEJ- Instead, we pick up a specific non-perturbative mechanism of the baryon 
number non-conservation, associated with triangle anomaly. The choice of this mechanism is 
explained, partially, by the authors' personal taste. In addition, anomalous fermion number 
non-conservation is a general phenomenon for theories with chiral fermions, and is present, e.g., 
in the standard model of electroweak interactions. This mechanism, being operative at high 
temperatures, may lead to the baryogenesis at the electroweak scale. 

The possibility that baryon asymmetry may be due to physics which is probed at accessi- 
ble energies has attracted a lot of attention recently and serves as a powerful motivation for 
the development of high temperature field theory, theory of phase transitions, non-equilibrium 
statistical mechanics. 

The fact that baryon number is rapidly violated at high temperatures |25[] (for earlir discussion 

see 



31| , D, |32|) and under other extremal conditions |33|, [34], [35], [36], |37|, |3^| naturally leads one 



to enquire whether electroweak baryon number non-conservations occurs at high enough rate in 
collisions of energetic particles. This problem has attracted considerable interest in recent years, 
after the first — and encouraging at the time — quantitative results were obtained ||, g|. In 



spite of remarkable theoretical developments, this problem is still not completely solved; existing 
results indicate that the electroweak baryon number violating processes occur at unobservable 
rates even at very high energies. 

The paper is organized as follows. In Section 2 we provide the necessary background and 
discuss the mechanism of anomalous non-conservation of fermionic quantum numbers together 
with relevant bosonic classical solutions (instantons and sphalerons). Section 3 contains prelim- 
inary discussion of the role of baryon number violating electroweak processes in early Universe. 
The fermion number non-conservation at high temperatures is considered in Section 4. In Sec- 
tion 5 we present recent developments in the theory of high temperature phase transitions. In 
Section 6 we briefly address the question of survival of the primordial baryon asymmetry. The 
discussion of various electroweak baryogenesis mechanisms is contained in Section 7. We turn 
to electroweak baryon number non-conservation in particle collisions in Section 8. Section 9 



4 Again, there are exotic mechanisms making use of the inflationary stage of the Universe expansion, in which 
the underlying theory conserves C or CP, the Universe as a whole is charge symmetric, but the visible part is 
not, see review [O. 
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contains concluding remarks. 



2 Basics of anomalous non-conservation of fermion quan- 
tum numbers 

Let us discuss non-perturbative non-conservation of fermion quantum numbers in the context 

of a model with the gauge group SU(2) and the massless left-handed fermionic doublets 

i = 1, . . . , til- The absence of global anomaly |41| requires that is even. We also add a Higgs 



doublet that breaks the SU(2) symmetry completely. Then this theory is a simplified version 
of the electroweak sector of the minimal standard model. All relevant features of the standard 
model are present in this simplified theory; later on we shall comment on minor complications 
due to U(l)y gauge symmetry, right-handed fermions and Yukawa interactions leading to fermion 
masses. One may regard the simplified theory as the standard model in the approximation where 
sin#iy and all fermion masses are set to zero; for three families of quarks and leptons one has 
n L = 12 and 

4 ] = ti f L a Jl}, (2-1) 

where / = 1,2,3 is the family index and a = 1,2,3 labels the colour of quarks. 
At the classical level, there exist til conserved global U(l) currents, 

which correspond to the conservation of the number of each fermionic species. At the quantum 
level these currents are no longer conserved due to the triangle anomaly f42, 43, 44 



<V£° = ^Tr (F, V F, V ) . (2.2) 

Therefore, one expects that fermion numbers Np = J d?x Jq are not conserved in any process 
where the gauge field evolves in such a way that 

N[A] = / d 4 x TtF, u F, u ^ 0. (2.3) 

Namely, 

AiVjP — N[A], i = l,...,n L . (2.4) 

It is clear from eq. (|2.3j ) that in weakly coupled theories, one has to deal with strong fields: 
the field F^ u = ^T a F^ v should be of order 1, and A® = O^g' 1 ). So, it is natural that the 
(semi) classical treatment of bosonic fields is often reliable. 

Equation ( |2.4j ) may be viewed as the selection rule: the number of fermions changes by the 
same amount for every species. In terms of the assignment ( |2.1| ) it implies, in particular, 

AN e = AN IM = AN T = N[A], 

AB = --3-3-N[A], (2.5) 
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where the factor 1/3 comes from the baryon number of a quark, while the factor 3 • 3 is due to 
colour and number of generations. So, the amounts of non- conservation of baryon and lepton 
numbers are related: 

AN e = AiV^ = AN T = ^AB, 

(B — L) is conserved while (B + L) is violated. 

The analysis of gauge field configurations with the non-zero topological number (|2.3| ) is con- 
veniently performed in the gauge 

A = 0. 

In this gauge, there exists a discrete set of classical vacua, i.e. pure gauge configurations 

<j) = LO<f) , 

where 0o = (0,t>/\/2) is the Higgs field in the trivial vacuum. The gauge functions u> depend 
only on spatial coordinates uj = w(x) and are characterized by an integer: 



/ d 3 x e^TrfudiU' 1 ■ wa-w" 1 • ud k uj- 1 ) 
2Att 2 J 



The vacua with different n[u] cannot be continuously deformed into each other without generat- 
ing non-vacuum gauge fields, so these vacua are separated by a potential barrier. Therefore, the 
gauge-Higgs system is similar to a particle in periodic potential, as shown in fig. [I]. An explicit 
construction of the minimum energy path connecting the neighboring vacua was carried out in 



ref. |[45|| , and the fermion sea contribution to this path was evaluated in ref. |46 
The topological number density entering eq. fl2.3p is a total derivative, 



2 Tr (yF^uF^u^ d^K^j 



32vr 

where 

K, = e^ A "Tr (f uX A p - ^A U A X A, 
If one is interested in vacuum-vacuum transitions, then 

t=+oo 



N[A]= J d 3 xdt d^K^ = J d 3 x K 



t=— oo 



= n[u t=+OQ ] - nfc^-oo]. 

So, the topological number of the gauge field is non-zero for transitions between the distinct 
vacua. 

At zero energies and temperatures, the transition between vacua with different n[u] is a 
tunnelling event which is described by instantons (constrained instantons in theories with 
the Higgs mechanism |47[| ). In pure Yang-Mills theory an instanton is the solution to the 
Euclidean field equations which is an absolute minimum of the Euclidean action in the sector 
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N[A] = 1. Properties of instantons are reviewed in ref. [48]. The instanton field, up to gauge 



transformations, is 



1 2x u 

ft ~ n ^ va ^1 i n 2 ' 



g x z + p 



(2.6) 



where rj^ va are the 't Hooft symbols, and p is an arbitrary scale to be integrated over. The 
instanton action is 



S, 



inst 



57T 



and the tunnelling amplitude is proportional to 



A ■ r\~ p Sins 
-tiinst ^ c 



In the electroweak theory, the tunneling probability is unobservably small, 



a inst oc exp 



4-7T 



aw 



10 



-170 



(2.7) 



(2.8) 



where 



a 



1 

29' 



sin 9 W 

In theories with the Higgs mechanism there appears a slight complication. There are no 
solutions to Euclidean field equations, i.e. no exact minima of the Euclidean action in sectors 
with N[A] 7^ 0. The reason is that the action for configurations like Q2.6|) , with appropriate 
Higgs field, depends on the instanton size p and decreases as p tends to zero. To evaluate the 
functional integral in that case one introduces a constraint that fixes the size of the configuration 



47[[ , then minimizes the action under this constraint and finally integrates over p. The outcome 



of this procedure is as follows. The instanton contribution into the functional integral becomes 

(2.9) 



f 4 dp ( 8tt 2 
J d x — p{p) exp I 



2 2 2 
7T V p 



where xo is the instanton position and p(p) is a function of g and p that varies relatively slowly. 
The integral (|2.9|) is saturated at p < v^ 1 , so that the size of the constrained instantons is 
smaller than the inverse Vy-boson mass mw = gv/2. The constrained instanton configuration is 
conveniently described in the singular gauge where the original pure Yang-Mills instanton has 
the form 

A*-U 2p2xU 
» g V( * ua x 2 (x 2 + p 2 )' 

The constrained instanton is given by eq. (p . 1 0|) at x m w , 

x, 

A inst oc Q~ mw \ x \ 

Clearly, the tunnelling rate is still suppressed by the exponential factor 

In this paper we discuss processes at high temperatures or energies. The relevant energy 
scale is set by the height of the barrier between different vacua as sketched in fig. [I]. This 



(2.10) 

and exponentially decays at large 
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height is determined by the static saddle point solution to the Yang-Mills-Higgs equations, the 
sphaleron j49], |32|. This solution was found previously in refs. |5T], [52], |53|, but its relevance 
to topology was realized only in ref. []32]| . By simple scaling one obtains that the static energy 
of the sphaleron solution in our simplified model, which is equal to the height of the barrier at 
zero temperature, is 

E , ph = ^ B f!M.)_ (2 . U ) 

aw \mw / 

where m# is the mass of the Higgs boson. The function _B(m#/mw) has been evaluated numer- 
ically P2 ; it varies from 1.56 to 2.72 as m#/my^ varies from zero to infinity (]. So, the height of 
the barrier in the electroweak theory is of order 10 TeV. 

At energies above E sp h the system can in principle evolve from a neighbourhood of one vacuum 
to another in a classical way, without tunnelling [| as outlined above, this classical process will 
lead to non-conservation of baryon and lepton numbers. Clearly, having enough energy is a 
necessary, but not a sufficient condition for the absence of exponential suppression of the baryon 
and lepton number violation rates. Whether the exponential suppression actually disappears or 
not is a matter of complicated dynamics which is one of the main subjects of this paper. 

There are at least two possible ways to see that fermion quantum numbers are indeed violated 
in instanton-like processes. One of them [^, |21[ makes use of zero fermion modes in Euclidean 
background fields with N[A] ^ 0. This approach is reviewed in ref. 



481 and its Minkowskian 



counterpart is considered in ref. |59[ . A more intuitive way |60| , |59|j is related to the phenomenon 
of level crossing, which is as follows. Consider left-handed fermions in the background field 
A(x, t) which changes in time from one vacuum at t — — oo to another vacuum at t — +oo (we 
again use the gauge A — 0). At each time t one can evaluate the fermionic spectrum, i.e. the 
set of eigenvalues of the Dirac Hamiltonian in the static background A(x, t) where t is viewed 
as a parameter. The spectrum varies with t; some levels cross zero from below and some cross 
zero from above. The relevant quantity is the net change of the number of positive energy levels, 
which is the difference between the total number of levels that cross zero from above and from 
below in the course of the entire evolution from t = — oo to t = +oo. A general mathematical 
theorem pT| says that this difference is related to the topological number of the gauge field, 



A. 



A_ = n[u t=+OQ } - n^-oo] = N[A\. 



(2.12) 



Recall now that at vacuum values of A, the ground state of the fermionic system has all negative 
energy levels filled and all positive energy levels empty. A real fermion corresponds to filled 
positive energy level and antifermion is an unoccupied negative energy level. As the energy 
levels cross zero, the number of real fermions changes, and the net change in the fermion number 
of each left-handed doublet is 



AN 



« 



N: - A_ 



Combining this relation with eq. ( |2.12| ) we see that the fermion number is not conserved indeed, 
and the amount of non-conservation is in perfect agreement with the anomaly relation, eq. ([2 . 1|) . 



5 At very large tuh /rriw the situation is more complicated |34], , but the estimate ( |2.1l[ ) remains valid. 

6 Of course, this travel does not typically proceed exactly through the sphaleron configuration. Some (not 
necessarily small) deformations of sphalerons are considered in refs. |fi6| , |57| , |5§|| where topological properties of 
these configurations are investigated. 
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Although the above discussion was for massless fermions, the results remain valid for the 
standard electroweak theory where fermions acquire masses via the Yukawa coupling to the 
Higgs field |63|, |64| . Indeed, the triangle anomaly for baryon and lepton currents remains 
valid in the standard model, so that relation Q2.5p must hold. The counting of fermion zero 
modes in the instanton background confirms this expectation [|62], |64], |65|. Also, the level crossing 



phenomenon has been explicitly found in theories of this type B7, 55 1. So, the complications 
due to right-handed fermions and fermion masses do not change the picture of baryon and lepton 
number non-conservation. 

Finally, the presence of U(l)y gauge symmetry in the standard model does not modify the 
analysis to any considerable extent either. There are no instantons of the U(l)y gauge field, 
while the effect of the C/(l)y interactions on the measure /x(p) for SU{2) instantons in eq. (pl9|) 
is tiny. Also, the energy of the SU(2) sphaleron is still given by eq. fl2.11|) where the factor B 
depends also on sin 2 W . For the actual value sin 2 Q w = 0.23, the deviation of B from its SU (2) 



values is numerically small [32, 68, 69 



3 Baryon asymmetry: preliminaries 

In this section we qualitatively discuss the issues relevant to the main topic of this review - 
electroweak baryon number non-conservation at high temperatures and generation of the baryon 
asymmetry of the Universe. These issues will be considered in much more detail in the following 
sections, so this section may be regarded as a guide for a reader not familiar with the subject. 
Most of what is said in this section should not be taken too literally: we will somewhat oversim- 
plify the picture of the electroweak physics in the early Universe and hence will use fairly loose 
terms. 

In hot Big Bang cosmology, there is an epoch of particular interest from the point of view 
of the electroweak physics. This is the epoch of the electroweak phase transition, the relevant 



temperatures being of the order of a few hundred GeV [TO], [H], |72|, ff3|. Before the phase transition 
(high temperatures), the Higgs expectation value is zero, while after the phase transition the 
Higgs field develops a non-vanishing expectation value. The critical temperature T c depends 
on the parameters of the electroweak theory; in the Minimal Standard Model (MSM) the only 
grossly unknown parameter is the mass of the Higgs boson, m#. In extensions of the MSM, 
there are more parameters that determine T c . 

At sufficiently small mu in MSM, the phase transition is of the first order, while at large mu 
the exact nature of the phase transition is still not clear: it may be weakly first order, second 
order or smooth cross-over. It is important that the masses of W- and Z-bosons immediately 
after the phase transition, mw{T c ) and mz{T c ), are smaller than their zero temperature values; 
the precise behaviour of mw(T) and mz(T) again depends on the parameters of the model (on 
mil in the MSM). Generally speaking, the stronger the first order phase transition, the larger 
m w{T c ) and mz{T c ). The electroweak phase transition is considered in more detail in section 5. 

Let us now turn to the rate of the electroweak baryon number non-conservation at high 
temperatures. While at zero temperatures the B non-conservation comes from tunnelling and 
is unobservably small because of the tunnelling exponent, it may proceed at high temperatures 
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via thermal jumps over the barrier shown in fig. |T] ||25|| . At temperatures below the critical 
one, T < T c , the probability to find the system at the saddle point separating the topologically 
distinct vacua is still suppressed, but now by the Boltzmann factor, 

/ E sph (T)\ 

Tocexpl — J, (3.1) 

where 

E vk (T) = 2 ^1 B (™L 
aw \mw> 

is the free energy of the sphaleron. Once the system jumps up to the saddle point (i.e. once the 
sphaleron is thermally created), the system may roll down to the neighbouring vacuum, and the 
baryon and lepton numbers may be violated. Therefore, the factor (|3.1| ) is also the suppression 
factor for the rate of the electroweak baryon number non-conservation at T < T c . 

At T > T c , the exponential suppression of the baryon number non-conserving transitions is 
absent. The power-counting estimate of the rate per unit time per unit volume in the unbroken 



phase is then 74, 75 



T = const • (a\yT) , (3.2) 

where the constant is of order 1. The rate of the electroweak B non-conservation is considered 
in detail in section 4. 

The rates (|3.1|) and ( |3.2| ) are to be compared with the rate of expansion of the Universe, 

1 ( T \ nn 

— = const • — — T, 
tu \M Pl J 

where the constant is of order 10 _1 . Clearly, in the unbroken phase the B non-conservation 
rate is much higher than the expansion rate in a wide interval of temperatures, T c < T < 
0.1 • Mpia w ~ 10 12 GeV. Therefore, the electroweak B non-conserving reactions are fast at 
these temperatures. After the phase transition the situation is more subtle: the rate of B non- 
conservation exceeds the expansion rate if the phase transition is weakly first order (mw(T c ) 
is small) or second order or of the cross-over type; on the other hand, the rate of I?-violating 
processes is much lower than the expansion rate if the phase transition is strongly first order 
(mw{T c ) is large enough). The electroweak B non- conservation switches off immediately after 

E (T ) 

the phase transition if ^ > 45 (see sections 6 and 7) and operates after the phase transition 
in the opposite case. This inequality is not satisfied in the MSM (section 7) with m top = 175 
GeV and experimentally allowed Higgs mass m# > 65 GeV. So, the 5-violating reactions are 
fast after the phase transition in the MSM. 

In the extensions of the MSM, the properties of the phase transition are determined by more 
parameters than just the zero temperature Higgs boson mass. So, for some region of the param- 
eter space, the electroweak B non-conservation is negligible after the phase transition. 

Clearly, the above observations are directly relevant to the problem of the generation of the 
baryon asymmetry of the Universe whose quantitative measure is the dimensionless ratio of the 
baryon number density to entropy density, 

= — ■ 
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This quantity is almost constant during the expansion of the Universe at the stages when baryon 
number is conserved, and its present value is 



-11 



A B = (4 - 6) x 10 

Several possibilities to generate the baryon asymmetry are discussed in the literature, which 
differ by the characteristic temperature at which the asymmetry is produced. 

(i) Temperature of grand unification, T ~ 10 15 — 10 16 GeV. 

A viable possibility is that the observed baryon asymmetry is generated by baryon number 
violating interactions of grand unified theories. The effect of the electroweak processes is basically 
that (B + L), generated at grand unified temperatures, is washed out at some later time (recall 
that {B — L) is conserved by anomalous electroweak processes). The asymmetry may survive 
from the grand unification epoch only if a large (B—L) asymmetry is generated at T ~ 10 15 — 10 16 
GeV, and there are no strong lepton number violating interactions at intermediate temperatures, 
100 GeV < T < 10 12 GeV (otherwise all fermionic quantum numbers are violated at these 
temperatures, and the baryon asymmetry is washed out). The first requirement points to non- 
standard, (B — L) violating modes of proton decay, though this indication is not strong. We 
discuss in section 6 some issues related to this scenario of baryogenesis. 

(ii) Intermediate temperatures, 1 TeV <T<C 10 15 GeV. 

An interesting possibility is that there exist lepton number violating interactions at interme- 
diate scales, and these interactions generate a lepton asymmetry of the Universe at intermediate 
temperatures. Then this lepton asymmetry is partially reprocessed into baryon asymmetry by 
anomalous electroweak interactions |75[]. Possible manifestations of this scenario are Majorana 



neutrino masses (which actually may be helpful from the point of view of solar neutrino exper- 
iments, for a review see, e.g., ref. J771) and/or lepton number violating processes like [i — > ej, 
fi — > eee and /i-e conversion. A more detailed discussion of this possibility, together with the 
analysis of concrete models, can be found in refs. [|7|, |(], [81], |82| . 

Another mechanism able to generate the baryon asymmetry at intermediate temperatures 



83j] deals with coherent production of scalar fields carrying baryon number. At a later stage 
the "scalar" baryon number stored in scalar fields is transferred into an ordinary one. The 
most recent consideration of this interesting possibility in the framework of the sypersymmetric 
standard model can be found in ref. ||B4|| . 

(iii) Electroweak temperatures, T ~ (a few) x 100 GeV. 

The remaining possibility is that the observed baryon asymmetry is generated by anomalous 
electroweak interactions themselves. Since the Universe expands slowly during the electroweak 
epoch, a considerable departure from equilibrium (the third Sakharov condition) is possible 
only from the first order phase transition. Indeed, this transition, which proceeds through 
the nucleation, expansion and collisions of the bubbles of the new phase, is quite a violent 
phenomenon. The dynamical aspects of the first order phase transitions in the Universe are 
considered in section 5. 

A necessary condition for the electroweak baryogenesis is that the baryon asymmetry created 
during the electroweak phase transition should not be washed out after the phase transition 
completes. In other words, the rate of the electroweak 5-violating transitions has to be negligible 
immediately after the phase transition. As outlined above, and discussed in detail in section 
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7, the latter requirement is not fulfilled in the Minimal Standard Model, so the electroweak 
baryogenesis is only possible in extensions of the MSM. Extending the minimal model is useful in 
yet another respect: it generally provides extra sources of CP violation beyond the Kobayashi- 
Maskawa mechanism, so that the second Sakharov condition is satisfied more easily Q The 
phenomenological consequences of these extra sources of CP violation are electric dipole moments 
of neutron and electron [R^, |8fifl , whose values are expected, on the basis of the considerations of 
baryogenesis j87|, to be close to existing experimental limits. 

Several specific mechanisms of electroweak baryogenesis are outlined in section 7. The outcome 
is that the observed baryon asymmetry may naturally be explained within extended versions of 
the Standard Model. This result is particularly fascinating as the physics involved will be probed 
at LEP-II and the LHC relatively soon. Naturally, most of our review is devoted to the topics 
related to the electroweak baryogenesis. 



4 Sphaleron rate at finite temperatures 

In this section we attempt to describe the present situation with the computation of the rate 
of fermion number non-conservation at high temperatures. We shall try to separate the exact 
results from (plausible) assumptions. We begin with the qualitative discussion of the rate and 
derive the Van't Hoff-Arrhenius type formulae for the rate, valid at sufficiently low temperatures. 
Then we derive an exact real-time Green function representation for the rate and show how it 
can be related to the more qualitative discussion. The quantum corrections to the rate are 
discussed. At the end, we present some numerical results for the sphaleron rate. 

Of course, there is much similarity between the description of sphaleron processes and reaction- 
rate theory in condenced matter physics. The latter is reviewed, e.g., in ref. 



4.1 Qualitative discussion 

As outlined in Section 2, the anomalous fermion number non-conservation is associated with 
the transitions of the bosonic fields from the classical vacuum of fig. [I] to the topologically 
distinct one. For the case of zero temperatures, low fermion densities and low energies of 
colliding particles, the initial state of the system as well as the final state are close to the vacuum 
configurations. So, to experience fermion number non-conservation, the system has to tunnel 
through the barrier. This process can be described by instantons and is strongly suppressed by 
the semiclassical exponent exp(— ^). 

In order to deal with topological transitions at non-zero temperatures let us consider first a 
simple example of the system with one particle in the double-well potential with the Lagrangian: 

L = h 2 -U(x), (4.1) 
U(x) = \\{x 2 - c 2 f. (4.2) 



7 Though it is still not excluded that the KM mechanism alone is sufficient for baryogenesis, see section 7. 
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The corresponding Hamiltonian is 



H = ip 2 + U(x), p = x (4.3) 



and the curvature of the potential at its minimun, x = c, is 

U"(x = c)=m 2 = 2 Ac 2 

Suppose that a particle is initially in the left well, and we want to calculate the probability 
of finding this particle in the other well. Let us take first the case of zero temperatures and 
consider the transition from the classical ground state. The probability of tunnelling can be 
found in the WKB approximation and is of the order of 

P ~ exp(-2S ), S = J C ^2U(x) dx. (4.4) 

It is exponentially suppressed provided that the energy barrier separating different classical 
ground states is sufficiently high. 

At finite temperatures, in addition to the ground state in the left well there are excited 
states with non-zero energy E. The probability of the state with the energy E is given by the 
Boltzmann distribution, exp(—E/T). Hence, the rate of the transitions is proportional to the 
sum of the probabilities of transitions from the excited levels with energy E weighted with the 
thermal distribution, 

P ~Y. e M-E n /T - 2S(E n )), (4.5) 

n 

where 

ME) , 

S(E)= dxJ2(U(x)-E), U{x{E)) = E. (4.6) 

J-x(E) V 

At temperatures T ^> m the sum can be approximated by the integral over x, 



dxexp 



_ 2 £ dyy/2(U(y)-U(x)) 



(4.7) 



with the result 



P ~ eM-Uo/T), (4-8) 

where Uq = U(0) = |Ac 4 is the height of the barrier. This result is clear from the physical point 
of view. Namely, P in eq. (4.8) counts the number of states with energy higher than the height 



of the barrier. At temperatures T < m the number of these states is exponentially suppressed 
by the Boltzmann exponent and their contribution is smaller than the contribution of tunnelling 
from the vacuum state. On the other hand, at high temperatures T ^> m the main contribution 
to the transition rate comes from the states with energy higher than the height of the barrier, 
which can overcome the barrier classically. Hence, we can address the problem of interest by the 
entirely classical calculation of the rate, which is equal to the probability flux in one direction 
(from left to right) at the point x = (see ref. |38j and references therein), 

/r/ w\-\ / dpdxexp(-H/T)5(x)6(p)p m . TT . . 

7 = {6(x)6(x)x) = F ; v ' = — exp {-U T ), (4.9) 

J dpdx exp(—H/ 1 ) 2tx 
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if T Uq. Note that the curvature of the potential near the saddle point at x = does not enter 
the final result; the quantum constant H does not appear in the answer at all. Note also that 
the classical treatment of the problem is applicable only if Uq/T < 2Sq. In the opposite case the 
rate of the quantum tunnelling is higher than the rate of the classical transitions. At the same 
time, the saddle point approximation we used for the calculation of the integral fl4.9|) is valid 
only for Uq/T ^> 1. If the latter relation does not hold, the calculation should go beyond the 
saddle-point approximation. This discussion can be easily generalized to the case of the systems 
with many degrees of freedom, in particular to the field theory we are interested in |89|, |X| . 

Let us consider specifically sphaleron transitions. As in the quantum-mechanical example 
discussed previously, we would like to put our system initially in the vicinity of one of the topo- 
logical vacua, say with n = 0, and determine the rate at which the system moves to neighbouring 
vacuum sectors. The sphaleron configuration, lying on a minimal energy path connecting two 
close-by vacua with different topological numbers, plays a crucial role in the computation. 

The energy functional near the sphaleron configuration can be written in the quadratic ap- 
proximation as follows, 

H = E sph - \ulxl + \ 5>?s» + i £p?, (4.10) 

where Xi and pi are the normal coordinates and momenta, Ui > are corresponding frequencies, 
and the index "-" refers to the negative mode. Now, the surface x_ = in the configuration 
space is the complete analogue of the saddle point x = in the 1 degree of freedom model 
considered above. If we put our system on this surface, and let it evolve with time, then it will 
almost definitely move to the sector n = 1 (0) (and stay there for long time) if the projection of 
initial momenta to the normal to the surface is positive (negative). So, to count the number of 
transitions with the topological number change we should calculate the probability flux through 
the surface X- = in one particular direction (cf. refs. [^, [K|). The rate per unite time and 
unite volume is given by WA pi], [75| 



1 1 r 

r = — — J dxidpidx-dp-5(x-)6(x~)x^ exp(— H/T), (4-H) 



where Z is the statistical sum 



Z = JvVVQexp(-H/T) (4.12) 



and V is the volume of the system. As in the simple example, the rate does not depend on 
the curvature along the sphaleron negative mode. After Gaussian integration over momenta the 
result may be written in a compact form 

r = Vzh^ F -"" (4 ' 13) 

where lmF sp h is the result of the formal computation of the imarinary part of the free energy 
near the sphaleron in the one-loop approximation (since the sphaleron is a saddle point of the 
energy functional rather than its minimum, the functional integral around it does not exist). 
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The treatment of the sphaleron zero modes is fairly standard. The total number of sphaleron 
zero modes is 6. Three translational zero modes restore the correct volume dependence of the 
rate while SU(2) transformations introduce some normalization factor. The final result for the 
rate reads WM 

Here the factors N tr — 26, N rot ~ 5.3 • 10 3 come from the zero mode normalization |74 |, k is 
the determinant of non-zero modes near the sphaleron. Again this result is purely classical, 
and it does not contain h. It is not applicable at very low temperatures where > S inst 
(there quantum tunnelling is more important than classically allowed transitions), and at high 
temperatures where the exponent is not large compared to 1. 

There are several important assumptions used in the above derivation of the sphaleron rate. 
The first one is the applicability of the classical theory to the description of the topology change 
in high temperature plasma. The second one is inherent in the classical theory itself. The 
thermodynamics of the classical field theory is, strictly speaking, ill defined due to the Rayleigh- 
Jeans instability (in field theory language there are ultraviolet divergences). Even besides this, 
we assumed that if the system is initially on the surface x_ = then it will move to the sector 
with topological charge n = 1 provided that p_ > and to n = in the opposite case Q. None of 
these asumptions has been proved with any rigour up to now. A number of quantum corrections 
to the rate, associated with the contributions of high momentum particles (k ~ T) can be taken 
into account, the infinities in the classical theory can be properly dealt with at the one-loop 
level, but the regular procedure of evaluating higher order (two-loop, etc.) contributions to the 
rate is not known. Moreover, the approach discussed above does not shed any light on the rate 
at high temperatures, when the sphaleron approximation breaks down. 

Below we will describe the Green's function approach to the problem of B non-conservation, 
which follows from the first principles of statistical mechanics and allows (at least in principle) 
the discussion of the B non-conservation rate beyond the semiclassical approximation j75], 0, |93| . 

4.2 The Green's function approach 

Let us consider first the behaviour of the quantity Q{t) 

Q(t) = f t q{x)d 4 x (4.15) 
Jo 

at high temperatures in a quantum system without fermions. Here q(x) is the topological number 
density 

If the system starts in the vicinity of one classical vacuum (say, n — 0), then, due to the 
sphaleron-like transitions, it will move randomly to the vicinity of the other vacua of fig. IB. 

8 This is only true in the quadratic approximation of the energy functional near the surface. Higher order 
terms in the expansion of the energy functional will, in general, introduce recrossings of the surface [ [T4J and 
modification of the rate. 
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Because of the periodicity of the static energy, these processes will not increase the free energy. 
In other words, the quantity Q makes a "random walk" in the space of configurations, and 



(Q 2 ) = 2TVt, t -> oo, (4.17) 

where V is the volume of the system, (O) = Tr(Op), and p is the equilibrium density matrix, 
p = 4 exp(— y)- The quantity T is nothing but the rate of the transitions with a change of the 
topological number per unit time per unit volume. It is given by the correlation function 



- lim lim / ((q(x)q(0) + q(0)q(x)))dTx. (4.18) 

2 t^oo V^oo Jo 



Now, we can derive a fluctuation-dissipation theorem showing that the kinetic coefficient 
describing the relaxation of the fermion number is directly related to the rate T. Let us add left- 
handed (and right-handed) massless fermions to our system, switch off the Yukawa interaction 
with scalar fields, and suppose that the only deviation from the thermal equilibrium is that 
assosiated with the presence of small lepton and baryon numbers. In other words, all other 
interactions are supposed to be faster than those associated with anomalous B non-conservation^. 
For simplicity we assume that the averages of all conserved fermion numbers are equal to zero 
0, so that the only non-vanishing global charge is ((B + L)) = nV. Here B and L are left 
baryon and lepton numbers (in the massless limit only left-handed particles participate in the 
anomalous processes). 

In principle, the time dependence of the baryon number can be found from the solution of the 
Liouville equation for the density matrix p(t), 

a -f=*lP,H] (4.19) 
with the following initial condition at t — 0: 



1 



po = — exp 



~{H + iio{B + L)) 



(4.20) 

where po is the initial value of the chemical potential and Z in the statistical sum. Then, 

(B + L)(t) = Tr[p(t)(B + L)\ = Tr[p (B(t) + L(t))\, (4.21) 

where B, L (B(t), L(t)) are the operators of baryon and lepton numbers in Schrodinger (Heisen- 
berg) representation. We expect that if the density n is small compared with T 3 , then the time 
dependence of n follows from the kinetic equation 

OTL 

= -r B n, (4.22) 

9 In reality this is not always the case. For example, the chirality breaking reactions for massive light quarks 
and leptons due to Yukawa couplings may be slower that the anomalous reactions. The modification of the 



kinetic equations in this case is discussed in refs. H, [95[ ^ §7). 
10 A more general case is considered in refs. [75|, 
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where is the rate of the fermion number non-conservation we are interested in. Probably, the 
easiest way to derive this kinetic equation is to make use of the Zubarev formalism of the non- 
equilibrium density matrix Zubarev defines two density matrices to be considered: the 
first one is the so-called local equilibrium density matrix, which is time dependent only through 
p(t), the instant chemical potential for the operator (B + L): 

p le = i exp + p(t)(B + L)]j . (4.23) 

The magnitude of p slowly varies in time due to B and L non-conservation. The average value 
of n is related to p(t) as follows, 

n = ±Tr[(B + L)p le ] = 2 -p{t)T 2 N f) (4.24) 

where Nf is the number of fermion generations and we take into account that the baryon number 
of a quark is |, and the number of colours is 3. The change in the baryon number is related to 
the change in the chemical potential, dn = ^NfT 2 dp(t). 
At the same time, due to the anomaly equation, 

^n(t) = 2N f Tr(q(t)p ). (4.25) 

The main difficulty is to compute the large time asymptotics of this expression. Zubarev argues 
that at t > Tb the following so-called non-equilibrium density operator can be used instead of 
Po- 



1 

Pzub = -x exp 



-j, (H + e J dt'e-^n(t')(B(t') + L(t')) 



(4.26) 



where operators B and L are taken in the Heisenberg representation and the limit e — > +0 is 
assumed. This density matrix is static in this limit. Now, expanding the density matrix (4.26) 
with respect to small /i, integrating over t' by parts and neglecting the time derivative of p one 
obtains 

Tr[q(t)p Zub ] = l -T 2 p(t)T B , (4.27) 
where the rate of the baryon number dilution is written in terms of the retarded Green function, 



* = ™ lir ?n / dt'd z x[q{t), n(t')]_ exp(-e(t - 0)- (4-28) 



Now, the use of the spectral decomposition for the correlation function ( 4.18|) and Green's 
function (|4.28 ) p2 , [93|| shows that these two functions are in fact equal to each other up to a 
coefficient containing the number of fermion generations. Finally, one finds the desired relation 

T B = 12iV/^p (4.29) 

If the Yukawa interactions with the Higgs particles are faster than the sphaleron transitions, 
then one obtains the coefficient -j- instead of 12 |95| . 
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4.3 The relation to the "probability flux" formulae 



At first sight, the correlation function describing the rate of topological transitions ([4.28D has 



nothing to do with the probability flux through the surface x_ = we found in the subsection 4.1. 
Below we will see that they are in fact the same in the Gaussian approximation to the classical 
theory |[75| , at temperatures below the sphaleron mass but above the particle masses. One may 
expect that the classical approximation to the correlation function may be good enough, since 
the quantum bosonic distribution functions at low momenta are the same as the classical ones. 



Consider for simplicity the purely bosonic theory. Notice that expression ( 4.18 ) allows for the 



naive classical limit h — > (we leave aside the question of renormalization for a moment), 

H(V, Q)' 



T dass = lim / VVVQexp 



T 



Q(t)q(0), (4.30) 



where H(V, Q) is the classical Hamiltonian depending on the generalized momenta and coordi- 
nates, q(0) is the density of the topological charge at time t — expressed through canonical 
coordinates and momenta V, Q, and Q(t) is the topological charge (^4.15|) where A(x,t) is the 
solution to the classical equations of motion with initial conditions V, Q. Derivation of the clas- 
sical limit goes precisely along the lines of the corresponding analysis for the quantum mechanics 
given in ref. ||100|| . 



Intuitively, the main contribution to the path integral ( |4.30|) comes from configurations that 
start in the vicinity of one of the classical vacua of our system, evolve in time, pass near the 
saddle point (sphaleron), and relax in the vicinity of the other vacuum. For these configurations, 
Q(oo) = ±1, depending on the direction of motion. 

Consider the classical configurations crossing the surface x_ = near the sphaleron at some 
time t\. They have the form: 

A(x) = A^x - x , fi) + ^2x n A n (x- x ,O), (4.31) 



'sph 



X - X , Q) + Xn<f>n(x ~ X , ft). 



Here xo and Q are the collective coordinates corresponding to the sphaleron translational and 
rotational zero modes, and Xi are small. According to the Liouville theorem, the phase space is 
invariant, so that we can write 

(VVVQ) (0) = (VVVQ) (tt) = n Vx n Vp n dxodnN tr N rat \p„\dti (4.32) 

where p n are the momenta corresponding to the coordinates x n , dx- = \p-\dti, N tT and N rot are 
the normalization factors for the translational and rotational zero modes. Now, 

J d 3 x J*' dhq(0) = J d 3 x J * dtig(-x , -t x ) (4.33) 

is nothing but the topological charge of the configuration which at time t\ belongs to the surface 
x_ = 0. Then, if the momentum corresponding to the negative mode is positive, the configuration 
under consideration was evolving in time with the increase of the coordinate x_, producing on 
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average the topological charge |. In the opposite case the average topological charge is — |. So, 
for these configurations 

J d 3 x J 1 dtig(-x , = -sign(p_). (4.34) 
and Q(oo) = sign(p_). Finally, the correlation function ( |4.30|) is 

11 r 

r = -—N tr N rot J Vx n T>p n \p_\dp- exp[-H(x n ,p n )] (4.35) 

coinciding with what was found previously. One can see that the set of assumptions used in 
the computation of this correlation function is precisely the same as that for the estimate of the 
probability flux. 

4.4 Quantum versus classical rate 

In spite of the fact that we were able to write an exact quantum real-time correlation function 
describing the sphaleron rate, there are no regular methods which allow for the actual computa- 
tion. The Euclidean (Matsubara) field theory perturbative methods are of little help here, since 
the analytical continuation to the real time is necessary, which is hardly feasible by making use 
of the perturbation theory. 

The following arguments suggest that the leading quantum effects may be absorbed into the 
parameters of the classical theory, at least at sufficiently high temperatures. In the consider- 
ation of the sphaleron rate the most important role was played by the properties of the static 
configurations of the gauge and Higgs fields. They are certainly influenced by the precence of 
the high-temperature excitations. Nevertheless, it can be shown (see the discussion in Section 
5) that all static quantum temperature Green's functions for bosonic fields 0, Ai coincide, up to 
0(g 3 ) terms, with the static temperature Green's functions for the classical bosonic theory with 
the Hamiltonian 



H = Jd 3 x [^(x) 2 + P*(x)P(x) + \g%G% + (A$ ) t (A$) + U(<f>) 



(4.36) 



U(4>) = m 2 (T)&® + A(T)($ t $) 2 . 

Here E t and P are the momenta conjugate to the fields Ai and 0. The coupling constants g(T), 
A(T) and the mass m 2 (T) of the classical field theory can be found by well-defined perturbative 
prescription, to be discussed in more detail in section 5. The classical theory does not contain 
fermions, which are integrated out by the procedure of dimensional reduction (see section 5). 
Static classical Green's functions are finite, provided known one- and two-loop counterterms are 
added to the Hamiltonian ( [4.361) . The perturbative transition to the classical theory is possible 
in weakly coupled theories only. Moreover, the static correlation lengths in the theory must 



n To be more precise, the equivalence holds for any form of the kinetic part of the classical Hamiltonian, pro- 
vided it contains momenta only. Then for static Green's functions the integrations over momenta and coordinates 
are factorized. This ensures the equivalence of the classical high temperature equilibrium statistics and quantum 
3d zero temperature euclidean theory. 
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be large enough (/ ^> (7rT) _1 , I 3> (gT)^ 1 ) for the approximation to be valid. The Rayleigh- 
Jeans instability is nothing but the infinite renormalization of the vacuum energy in the 3d 
theory, which can be removed by vacuum counterterms. Therefore, we see that the static energy 
barriers can be found from the analysis of the saddle points of the energy functional defined 
in eq. ( 4.36 ). If m 2 (T) is negative, the symmetry is spontaneously broken, and the sphaleron 



solution does exist. Its energy now is temperature dependent via quantum corrections to the 
parameters of the classical theory. The size of the high-temperature sphaleron is of the order of 
the static correlation length I and is much larger, according to our assumption, than the inverse 
temperature (typical distance between particles). Now we recall that in the classical computation 
of the probability flux through the surface x_ = the integration over momenta is Gaussian (see 
eq.( 4.11|) ) and the main contribution comes from momenta ~ T. Therefore, for the normal 



sphaleron modes with uj n <C T the real time motion indeed can be considered as the classical 
one. This is not true for u n > T, but high energy sphaleron modes with u n ^> / _1 are close to 
those around the vacuum configuration and, therefore, cancel out in the rate computation. 

In the symmetric phase, at m 2 (T) > 0, the sphaleron solution does not exist. However, the 
typical static correlation length in the system, / ~ (g 2 T) _1 , is still much larger than the inverse 
temperature. It is natural to assume that the typical size of the configurations contributing to 
the rate is of the order of the correlation length; then the argument given above indicates the 
possibility of the classical description. 

So, the conjecture (which has never been proved, though) is that the rate of the fermion number 
non-conservation in quantum field theory at high temperatures (T ^> m^(T), gT 3> my/{T)) 
is given, up to terms 0(g 3 ), by the classical correlation function ( |4.18| ) with the Hamiltonian 



( [4.36| ). Since the statistical mechanics of the classical field theory does not exist due to ultraviolet 
divergences, the latter statement requires clarification. To define the classical statistics, an 
introduction of some high energy cutoff A is necessary, together with 3d counterterms removing 
divergencies from the potential part of the classical Hamiltonian. For this conjecture to be true, 
the necessary condition is that the classical correlation function ( |4.18| ) exists in the limit 

lim lim lim G (4.37) 

A-^oo t^co V-^oo 

(the order of limits is essential here, see below). Some evidence in favour of the last statement has 
been found by a number of direct real time Monte Carlo simulations, which will be considered 
in the last subsection. If the conjecture is true, it allows an immediate determination of the 
parametric dependence of the rate of sphaleron transitions: 

r = ( OW T)V(p,^), (4-38) 

since the rate has dimensionality (GeV) 4 , and the classical dynamics of the theory is governed 
by the unique dimensionful coupling g 2 T and two dimensionless ratios. Deeply in the symmetric 
phase, at T ^> T c , the scalar degrees of freedom decouple, since m 2 (T) ^> g A T 2 . Then Y ~ 

(« w t) 4 [bih- 



Recently, another conjecture was put forward in ref. |101|]. The authors suggest that instead of 



the classical Lagrangian ( |4.36|) one should use the Braaten-Pisarski effective action ||102|| , which 
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sums up so-called hard thermal loops in the particle amplitudes. This action is a generating 



functional for the gauge and matter fields at soft momenta; it was used in ref. [ 103 1 for the 
calculation of the the damping rate of fermions in the plasma. Yet another suggestion is to use 
Langevin-type equations with a friction term together with a random force |]104j| , instead of any 



type of deterministic equation of motion. The random force is served to mimic the interaction of 
the classical soft momentum modes with the short wave quantum fluctuations. These conjectures 
remain not proved either. 

4.5 The sphaleron rate in the broken phase 

The discussion in the previous subsections shows that in the broken phase, the rate of sphaleron 
transitions is given by expression ( |4.35| ), where the classical static energy which should be used 
in the evaluation of the functional integral is given by eq. ( |4.36| ), with temperature-dependent 
masses and couplings. So, in this regime, the rate is given by eq . ( |4. 14j ) with the replacement 
M w — > M W (T) where M W (T) = ^g(T)v(T), and the expectation value of the Higgs field is to 
be determined by the minimization of the classical potential U((f)). 

To completely determine the rate, one should calculate the 3-dimensional determinant of small 
fluctuations around the sphaleron solution, which is also temperature-dependent. Formally, it 
diverges linearly with ultraviolet cutoff, but this divergence is removed by the one-loop countert- 
erm for the scalar mass m 2 (T). Recently, this computation has been done numerically in refs. 



109| , |110|, |1 1 1|] where the values of the determinant can be found at different scalar coupling 
constants. For the small ratio of X(T)/ g 2 (T) the result has a simple form. Namely, instead of 
taking the tree value for the vacuum expectation value for the scalar field, one may obtain it 
from the minimization of the one-loop effective potential, 

1 / 6m 3 T + mf + 3ml) , (4.39) 



12tt 

where the mean field-dependent masses are defined as 

m T = ^g{T)<f>, m\ = m 2 M + 3A(T)0 2 , m\ = m 2 ^) + A(T)0 2 . (4.40) 



Then the 3d determinant k ~ 1; the precise numbers are given in ref. ||111|| . This is the most 
complete calculation of the sphaleron rate in the framework of 3d approach done up to now. 

Recently both the bosonic and fermionic determinants in the background of the sphaleron have 
been calculated for any temperatures in refs. ||105| , |106|1 . Probably, this is the most involved and 



non-trivial computation done until now for the sphaleron rate. Also, the change of the fermion 
number and behaviour of the relevant fermion level during the sphaleron transition has been 
followed. The authors conclude that the fermionic contribution is numerically important for 
m t = 175GeV. However, those calculations become unreliable close to T c [^} Fortunately, the 

12 We thank D. Diakonov and K. Goeke for discussion of this point. 
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full theory in the vicinity of T c can be reduced to a 3-dimensional problem, and one can show 
that almost entire effect of the determinants can be absorbed into the definition of coupling 
constants of the effective 3-dimensional theory (see section 4.4 and ref. [ 1Q7| ). If one uses these 



couplings to redefine the sphaleron solution, the remaining contributions are shown to be small 



108|| . As a result, the sphaleron rate in the Higgs phase is reliably calculable at the critical 
temperature for the theory with the scalar self-coupling up to \{T)/g 2 {T) < 0.04, see discussion 
in Section 7.1. 

The evaluation of further corrections to the rate is a problem which has not been solved, even 
the strategy of the necessary computation is not known. Parametrically, they are 0[g 2 T/Mw{T)}. 

4.6 Real time numerical simulations 

Making use of the conjecture on the possibility to calculate the quantum sphaleron rate within 



the classical field theory, one can compute the rate by the numerical simulations [|112j , |113| , [I I I 
The discretization of space, necessary for the numerical methods, provides a natural ultraviolet 
cutoff. Probably, the most convenient discrete formulation is provided by the lattice gauge 
theories. 

The Monte Carlo (MC) numerical computation of the correlation function (Q 2 (t)) consists 
of two steps. First, one should generate a set of configurations (coordinates and momenta) in 
accordance with the Boltzmann distribution exp(— if/T). Then, these configurations are used as 
initial conditions for the classical equations of motion. These equations are solved numerically, 
and the topological charge is computed as a function of time. Finally, the averaging of the 
quantity Q 2 (t) is performed. The computation should be repeated for different volumes of the 
system and for different lattice spacings; the extrapolation to the continuum limit is to be 
performed at the end. Many details of the described procedure can be found in the original 
papers JTTJ, [TJJ, [TT5], [TTB], |TT7], [TTB|, |T2{], [T2TJ, here we just present some results. 



Up to now, real time MC simulations of the topological number change were performed in 
three different gauge field theories. The first one is the 1 + 1 dimensional U(l) Higgs model, the 
second is SU(2) gauge-Higgs system, and the third is a pure Yang-Mills SU(2) theory. 

£7(1) theory in 1+1 dimensions. The Lagrangian of the U(l) Higgs model has the form: 



L = -\F^ V + (£> M 0)t(£> M 0) - V(<j>), (4.41) 



where F^ u = d^A v —d u A^ is the gauge field strength, is the charged scalar field, = d^—ieA^, 
and V((f>) is the scalar potential, 

K(0) = ^-c 2 ) 2 . (4.42) 
The topological number in this theory in (1 + 1) dimensions is 

Q(t) = ^J cfxt^F^. (4.43) 
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The static high temperature effective one- dimensional action is just [(Di<f>Y(Di(f)) + V(0)]. It 
corresponds to finite field theory in one dimension. The classical statistical mechanics of this 
system is described by the Hamiltonian 



-E 2 + PP* + (D 1 <py(D l (f>) + V((j>) 



(4.44) 



together with the Gauss constraint 



d x E - ie(0*P - 0P*) = 



(4.45) 



imposed on the admissible states. 

Let us put this system in a one-dimensional box of length L, —L/2 < x < L/2 and impose 
periodic boundary conditions on the fields. Then in the limit L — > oo the saddle point of 
this action - sphaleron - is a gauge transformation of the usual kink of the scalar field theory 
TT2| , |TT3| PI : 

(4.46) 
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The high temperature sphaleron rate in this theory has been calculated in refs. [91| |122| 
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(4.47) 



(4.48) 



where the sphaleron mass is given by E sp h = :y y^c 3 . The function f(x) for large x has been 
evaluated in ref. |91| (f(x) = ^-2 2x ~'i) and, for arbitrary x, in ref. ||122j| . 

frn. 



The real time dynamics of the classical system was studied on the lattice in refs. 



120, 121 



The results of the real time MC simulations show that the rate of sphaleron transi- 
tions, indeed, does not depend on the lattice spacing, provided it is small enough. Moreover, 
quantitative agreement with one- loop formulae is found at T < E sp h. At very high temperatures, 
T > E sp h the rate cannot be calculated analytically, and only the numerical results exist there 
123 ]. The numerical simulations with Langevin-type equations replacing the real time dynamics 
can be found in ref. ||104|. 



SU(2) gauge— Higgs system. In refs. ||1 1 5| , |116|| the study of the sphaleron transitions in 
the symmetric phase of the SU(2) gauge-Higgs theory was performed. The sphaleron transitions 
were clearly observed with different lattice spacings and volumes. However, the quality of the 
lattice data did not allow the extrapolation of the results to the continuum limit. All the data 
are consistent with the rate T = k(«vi/T) 4 with k > 0.4 0. 

SU(2) pure gauge theory. Recently, the systematic study of the sphaleron transitions has 
been performed in a pure SU(2) theory |119| . The use of the pure Yang-Mills theory instead of 
the full gauge-Higgs system is legitimate at the temperatures far above the critical one. At these 



13 In ref. [116] the corresponding constraint reads as k > 0.1; an arithmetic error of a factor 4.4 in this estimate 
was corrected in ref. [124]. 
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temperatires the scalar fields have masses ~ gT and decouple from the long-range fluctuations 
k ~ g 2 T which are believed to govern the topology change. The 3d part of the classical model is 
a field theory free from ultraviolet divergences. The authors show nice evidence that the classical 
rate exists (i.e. that it does not depend on the volume, if it is large, and lattice spacing, if it is 
small). Numerically, 

T = K(a w T) 4 J k = 1.09 ± 0.04. (4.49) 

Near the phase transition the scalar degrees of freedom have masses ~ g 2 T and do not decouple 
from the gauge fields. However, since the parameter g 2 T provides the only dimensionful scale of 
the problem, the parametric dependence of the rate is still the same, but the numerical coefficient 
may be different. 



4.7 Strong sphalerons 

We include in this section also the discussion of other high temperature processes, associated 
with anomaly, but now in QCD. They change the chirality of quarks and, since anomalous 
.B-violation deals with left-handed fermions, may be relevant for the discussion of the baryon 
asymmetry. 

It is well known that the quark axial vector current has an anomaly and therefore is not 
conserved. The rate of chirality non-conservation at high temperatures T str0 ng is related to the 
rate of topological transitions in QCD ("strong" sphalerons 125fl ), 

dQ 5 _ 72 

where Q$ is the axial charge. By analogy with the weak sphalerons, the rate of the strong 
sphaleron transitions is given by 

T strong = K SUi?) (a s TY, (4.51) 
where Ksu(3) is an unknown pure number. The characteristic time of these transitions is therefore 

Tstron9 = 72K SU(3) air (452) 

Taking as an estimate Ksu(3) ~ 1-3, this time is of order 30 ^! 00 , i.e. the rate of strong sphaleron 
transitions is comparable to or even higher than the rate of chirality-flip transitions mediated 
by the Yukawa coupling of the top quark. 



4.8 Concluding remarks 

It is by now well established that there is no suppression of the fermion number non- conservation 
at high temperatures. However, the quantitative formalism allowing for the calculation of the 
rate beyond the lowest order semiclassical approximation is still lacking. From our point of view, 
the most important challenge here is to establish the relation between the quantum rate and 
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the classical rate; the latter can then be computed with some kind of MC numerical analysis. 
Even in the framework of the classical physics it would be important to have an analytical 
understanding of the fmiteness of the rate in the continuum limit. Of course, the real time MC 
simulations in the broken phase of the SU(2) gauge-Higgs system would be very important, in 
particular because the rate is now known in the Gaussian approximation. 



5 Phase transitions in gauge theories 

Potentially, phase transitions provide a source of deviations from thermal equilibrium in the 
early Universe. Usually, in the theories with scalars (such as grand unified theories or the 
electroweak theory) the symmetry is restored at high temperatures and, at low temperatures, it 
is broken |70|, |7TJ . If the phase transition is of the first kind, it proceeds through the nucleation 
of bubbles of a new phase ||1 26| , |127|| . Depending on the parameters of a model, this process 



can be quite violent, the motion of the domain walls disturbs the plasma and may trigger the 
baryon asymmetry generation. In order to have the detailed non-equilibrium picture of the 
phase transition, a number of very difficult problems are to be solved. The questions potentially 
interesting for baryogenesis include the bubble nucleation rate, the structure of the domain walls, 
their velocity, the distribution densities of particles near the domain walls, etc. It is hard to 
get reliable answers to these questions, since they all deal with complicated non-equilibrium 
phenomena. Moreover, even the equilibrium treatment of the phase transitions faces a number 
of difficulties, associated with the so-called infrared problem in thermodynamics of the gauge 
fields. 

There are many excellent reviews and books devoted to the study of the phase transitions in 
gauge theories, see, for instance ref. ||128| , |129| , |130| , pL31f] . Our purpose in this section is to report 



on the progress that was achieved in this area in the last few years. Our special interest is in the 
study of the first order phase transitions which are strong enough to suppress the B-violating 
reactions in the Higgs phase (see Sections 3 and 7.1). The dynamics of much weaker phase 
transitions may be different and is a topic of extensive studies of refs. |pL32 , |133| , |134| , |135|| . 



In the first subsection we shortly discuss the validity of the equilibrium approximation to the 
description of the phase transition and present some useful equations for the determination of the 
bubble nucleation rate. Then we review the "rules of the game" , which allow an estimate of the 
relevant parameters of the phase transition with the help of a unique function - the perturbative 
effective potential for the Higgs field. Everything contained in these subsections has been known 
for a long time and is presented here for the sake of completeness. In subsection 5.3 we explain 
why the perturbation theory fails to describe the high temperature phase transitions. In the 
following subsections (5.4 and 5.5) we describe the formalism which allows us to determine 
reliably the parameters of the phase transition in weakly coupled theories. The specific results 
for the electroweak phase transition are discussed in subsection 5.6. Subsection 5.7 is devoted 
to the dynamics of the phase transition. 
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5.1 Equilibrium approximation 

Let us take for simplicity the Minimal Standard Model of electroweak interactions and consider 
it in the cosmological context. Suppose that the temperature of the system, T, is of the order 
of the VT-boson mass - the scale at which the electroweak phase transition is expected to take 
place. The first question that arises is whether the equilibrium description of this system is at 
all possible in the expanding Universe. In order to check this, we may compare the rates of 
different particle reactions in the standard model with the rate of the Universe expansion, t v , 
where the Universe age tjj is given by 

tu = (5.1) 

Here M = M P j/1.66iV* ~ 10 18 GeV, M Pl is the Planck mass, and N is the effective number 
of the massless degrees of freedom. The expansion rate of the Universe is the unique non- 
equilibrium parameter of the system before or some time after the phase transition; during the 
phase transition another typical non-equilibrium time scale, associated with the motion of the 
bubble walls, is relevant. This time scale is smaller than the Universe expansion rate by many 
orders of magnitude (see below) and thus the deviations from the thermal equilibrium are much 
stronger. 

Before or after the phase transition, the fastest perturbative reactions are those associated 
with strong interactions (e.g. qq — > GG); their rate is of the order of (T stron£ ,) _1 ~ (anv) ~ Oi 2 s T. 
Here a is the cross-section of the reaction, n is the particle concentration, v is the relative velocity 
of the colliding particles. The typical weak reactions, say ev — > eu, occur at the rate (T TOeo fc) _1 ~ 
aiyyT, and the slowest reactions are those involving chirality flips for the lightest fermions, e.g. 
e^H — > vW with the rate (T e ) _1 ~ f^awT, where f e is the electron Yukawa coupling constant. 
Now, the ratio -r- varies from 10~ 14 for the fastest reactions to 10~ 2 for the slowest ones; this 
means that particle distribution functions of quarks and gluons, intermediate vector bosons, 
Higgs particle and left-handed charged leptons and neutrino are equal to the equilibrium ones 
with an accuracy better than 10 -13 ; the largest deviation from thermal equilibrium (~ 10" 2 ) is 
being expected for the right-handed electron. 

These estimates show that the equilibrium description of the system is a very good approxi- 
mation before the phase transition, and soon after it is completed [5 Moreover, since the phase 
transition is expected to proceed through the bubble nucleation, the equilibrium description of 
the plasma is possible in the regions far enough from the moving domain walls. 

The above remarks suggest that the equilibrium statistical mechanics may be applied for the 
evaluation of the "static" properties of the phase transitions, such as the critical temperature T c , 
latent heat, jump of the order parameter (expectation value of the scalar field), interface tension 
(surface energy density of the plane domain wall separating different phases) . Another important 
characteristics is the metastability range: at the upper T + (lower T_ ) spinodial decomposition 
temperatures broken (symmetric) phase ceased to exist as a metastable state (fig. 0). The static 
correlation lengths for various operators may help to understand the structure of the domain 



14 Clearly, this is a model-dependent statement. For instance, if the Universe was as hot as, say, 10 17 GeV, 
then the equilibrium description of the (grand unified) phase transitions would be questionable, since the ratio 
7 s - would be of the order of 1. 

tu 
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walls. 

The first order phase transition in the early Universe is not an instant process, and its total 
duration is of the order of the Universe age. The bubbles of the new phase start to nucleate 
at temperatures somewhat lower than T c ; the bubbles expand and finally fill out the Universe 
with the new phase. This happens at temperatures above T_ . The fraction of the volume of the 
Universe occupied by a new phase, V(t), can be estimated in the following way [ |136j| . Suppose 
that the bubble nucleation rate per unit time per unit volume is 7Z(T), and the velocity of the 
bubble walls is constant]^] and equal to v. Then, 

P(t) = l-exp(-A(t)), (5.2) 

where V(t,t ) = ^-v 3 (t — t ) 3 is the volume that the bubble nucleated at time t Q occupies at 
time t: 

A(t) = fdt V(t,t )K(T(t )) (5.3) 
Jt c 

and t c is the time corresponding to the critical temperature, T(t c ) = T c . This equation does 
not take into account the red shift of the bubble velocity, which is in fact unimportant at the 
electroweak scale. Introducing the variable 

T c -T 

x = 



T 

± c 

and assuming that it is small (this is satisfied in the electroweak case), one obtains a simplified 
expression 

A=»^Vf^W. (5.4) 



3 \TJ Jo T c 4 



The phase transition is completed when 

A~l. (5.5) 

Since the electroweak scale is much smaller than the Plank scale, the probability of the bubble 
nucleation 1Z at the percolation time is very small: 

„K{T) (T c 



The typical bubble size is of the order of R ~ 2vjgx, where x is to be found from eq. ( |5~5l) . 

Computation of the bubble nucleation rate in the general case is a very complicated prob- 
lem, which has not been solved. The reliable estimates exist only in the so-called thin wall 
approximation, and the leading contribution can be read off from the Landau-Lifshitz book on 
statistical mechanics. Suppose that the temperature of the system is T < T c , and i < 1. Then 
the free energy of the critical bubble can be found from the minimization condition 

dF . . „ 4nR 3 . . , 

— = 0, F(R) = AnR 2 a - —Ap, (5.6) 



15 In the hot plasma, contrary to the vacuum case, there is a friction force acting on the bubble wall. This 
ensures the constant velocity of the wall. 
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where Ap = Lx is the pressure difference, L is the latent heat of the transition, and a is the 
surface tension. From the latter relation one immediately obtains 



K 



fcT c 4 exp( 



-S(x)), 



(5.7) 



where the bounce action is given by 



S(x) 



16na 3 1 
3T r L 2 1? 



[l + 0(x)) 



(5- 



and corrections 0(x) are model dependent. Indeed, the domain wall thickness is of the order 
of the typical correlation length £ in the system, which means that the radius of the bubble is 
defined up to corrections of order £. This produces an uncertainty in the action 8S/S ~ 
and gives an obvious requirement of the validity of the thin wall approximation R 3> £. The 
calculations of the bounce action in various models can be found in ref. ||131| , |137 . 

If R ~ £, then the thin wall approximation breaks down, the nucleation rate cannot be 
expressed only through macroscopic parameters of the phase transition (latent heat and surface 
tension) at the critical point. In this case the phase transition is delayed and the Universe is 
supercooled in the symmetric phase. The calculation of the bubble nucleation rate in this case 
cannot be done for a generic gauge theory because of an infrared problem in the thermodynamics 
of the gauge fields (see below), but it is feasible for weakly coupled pure scalar theory. A 
detailed study of the bubble nucleation in the scalar mean field theory defined by potential 
( p,13| ) is contained in ref. ||138|| . Naturally, the bubble nucleation rate receives dependence 



on the scalar correlation length £ at the phase transition, in addition to the surface tension 
and the latent heat. We refer here to [ |131| , |137| , |139| | for more details. The more complicated 
problem is to determine the prefactor in the expression for the bubble nucleation rate. Its 
computation in one-loop approximation in scalar models and electroweak theory has been done 
in refs. pO, MR p|, |143j. 



5.2 Simple estimates 

To describe the high temperature phase transitions in any given theory it is very important 
to have a relevant calculational formalism. The traditional tool is the effective potential for 
the scalar field 0. It is defined as the value of the free energy of the system (pressure with 
the minus sign) in a uniform background field 0. The minima of this potential correspond to 
the (meta)stable states of the system. The system undergoes a first order phase transition if 
there are two degenerate minima of this potential, separated by the energy barrier. In general, 
the effective potential is a gauge dependent quantity; perturbative calculations often produce 
complex terms. However, the values of the potential at its minima are gauge invariant; this 
allows for the gauge- invariant definition of the critical temperature and latent heat. 

The following simple strategy (the drawbacks of which we discuss later) gives a reasonable 
qualitative description of the phase transitions and often allows fairly accurate estimates [ |144|| : 

Step No. 1. Take your model and calculate the one-loop high temperature effective potential 
V{<f>,T). 
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Step No. 2. Define from it the critical temperature, jump of the order parameter, latent heat 
and surface tension with the use of the following equations. 
-T c and order parameter C : 

dV{ Q C ^ c) = 0> V(<f> c ,T c ) = V(0,T c ). (5.9) 
-Latent heat and surface tension: 



■ c 



(V{</> c ,T c )-V(0,T c )), (5.10) 



a = £ C y/2V{<j),T c ) #. (5.11) 

Step No. 3. Calculate the bubble nucleation rate in the thin wall approximation and compare 
it with the rate of the Universe expansion. Determine the bubble nucleation temperature and 
check the validity of the thin wall approximation. If it does not work, evaluate the bubble 
nucleation rate for a thick wall. To this end find 0(3) symmetric configurations extremizing the 
3d action 

(5.12) 



I(^0) 2 + \/(0,T) 



with the boundary condition <fi — > at x — > oo. The bubble nucleation rate is then 1Z ~ exp(— S). 

An example is provided by the Minimal Standard Model. Here the one-loop effective potential 
in the high temperature approximation is (for simplicity, we take the case when the Higgs boson 
is sufficiently light, and neglect the effects of the U(l)y interactions): 

V(<f>, T) = i 7 (T 2 - T 2 )0 2 - l -aT^ + ^A0 4 . (5.13) 
For the standard model with the top quark mass m top 

a = 9^/(3270,7 = ^ 2+ + A + _^ ; (5.14) 
and the lower metastability temperature T_ is related to the Higgs mass m 2 H = 2\v 2 through 

(5 - 15) 



Due to the presence of the cubic term, the potential predicts the first order transition with the 
critical temperature 

and the jump of the order parameter 
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Phase transition gets weaker when the scalar self-coupling increases. This is seen from the 
behaviour of the order parameter, latent heat, and the surface tension, all of which decrease 
with the increase of A. At large A the bubble nucleation rate can be determined in the thin 
wall approximation, while at small A (A ~ g 3 ) it breaks down, and the phase transition occurs 
with considerable supercooling. Qualitatively, the one-loop description gives correct results, but 
concrete numbers may be quite different from those obtained by a more refined treatment. The 
effect of the Debye screening on the effective potential was discussed in refs. ||145| , |137| , |146| , |147|| , 



and the two-loop computation has been done in refs. [|148|, |149| , |150| , |151| , |152| , |153| , |154| , |155|1 



Various aspects of the phase transition were discussed in refs. ||156| , |157| , |158| , |159 |. 



5.3 The infrared problem and factorization 



1601 , |T6T| that the per- 



It was realised by Linde and Gross, Pisarski and Yaffe a long time ago 
turbation theory for non-Abelian gauge theories with small coupling constants inevitably breaks 
down at high temperatures, at least in the symmetric phase. The physical reason is that at high 
temperatures, instead of the usual 4-dimensional expansion parameter g 2 the relevant parameter 
is p ~ g 2 riB(E), where ub(E) = [exp(—E/T) — is the Bose distribution function, E is the 
typical energy of a given process in the plasma. At E < T the expansion parameter is greater 

n 1rp 

than that at zero temperatures, namely, p ~ jf-> accounting for typical Bose amplification of 
the amplitudes. In the symmetric phase, gauge bosons are massless in perturbation theory, there 
is no infrared cutoff, and the expansion parameter can be arbitrarily large. In the broken phase 
the infrared cutoff is provided by the vector boson mass, and perturbation theory converges 
provided <C 1. Below we will give a more formal description of the infrared catastrophe. 

The fact that the perturbation theory breaks down at T ^ poses non-trivial difficulties for 
the description of the phase transition. Indeed, the phase transition occurs when the free energy 
of the broken phase is equal to that of the symmetric phase; but the latter cannot be calculated 
perturbatively. The latent heat of the transition receives contributions from both the symmetric 
and broken phases, and the same is true for the surface tension. In the Universe the phase 
transitions occur when it is cooling, so that the initial phase is the one in which perturbation 
theory breaks down. For strongly first order phase transitions, the leading contribution to the 
above parameters comes from the broken phase, where perturbation theory is applicable; in that 
case the perturbative description may be reliable. However, the infrared problem is fatal for 
an attempt of the perturbative quantitative study of the weakly first order phase transitions. 
Unfortunately, direct Monte Carlo lattice simulations of high temperature gauge theories are not 
possible at present for realistic theories, containing chiral fermions, due to well known difficulties 
in discretisation of the chiral fermion determinant. The purely bosonic sector of the models can 
be put on 4d lattice, and extensive 4d numerical simulations have been carried out in refs. 
|162| , |163| , |164j , |165| , |166|| , for a summary of results see a nice review by Jansen ||167|| . 



Recently, an approach was suggested, which allows for a solution of the equilibrium problem of 
phase transitions in weakly coupled (at zero temperatures) gauge theories ||152| , |168| , |107| , |169| . 
It combines both perturbative analysis and numerical Monte Carlo methods. The main idea 
of the method is the factorization of the different scales appearing in the description of the 
high temperature plasma. At the first stage, a much simpler effective theory, incorporating 
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all essential non-perturbative dynamics of the phase transition, is constructed by perturbative 
methods. The idea of this construction, known as dimensional reduction, goes back to the 
papers by Ginsparg JL7(J, and by Appelquist and Pisarski ||171| . It was developed in refs. 



152| , |168| , |107| , |169| in application to the phase transitions and applied later to hot QCD in 



ref. ||172| , |173| , |174| , |175| , |176|| . Different aspects of dimensional reduction were studied in refs. 



||177| , |178j , |179| , |180| , |181| , |182|| . At the second stage the effective 3-dimensional theory is analysed 



by non-perturbative methods (MC lattice simulations) 0. In the discussion below we follow ref. 



The idea of dimensional reduction comes from an observation that equilibrium finite temper- 
ature field theory is equivalent to Euclidean zero temperature field theory defined on a finite 
"time" interval j3 = 1/T supplied with periodic boundary conditions for bosons and antiperiodic 
ones for fermions. 

Periodic and antiperiodic boundary conditions enable one to decompose Bose (0) and Fermi 
(ifj) fields in Fourier series with respect to the finite time interval, 

oo 

<f ) (x,r)= J2 <Pn(x)exp(iuj b n T), (5.18) 



ip(x,r) = ^n(a;)exp(^r), (5.19) 

n=— oo 

where uo b n = 2nnT, u[ = (2n + 1)ttT. Therefore, 4d finite temperature field theory is equivalent 
to the 3d theory with an infinite number of fields, and 3d boson and fermion masses are just 
frequencies u b and . One can easily recognize here a perfect analogy to Kaluza-Klein theories 
with compact higher- dimensional space coordinates. The equilibrium dynamics of the theory 
is completely characterized by the set of Matsubara (imaginary time, or Euclidean) Green's 
functions, G n (u>i, ki), where Ui are discrete frequencies, n is the number of legs. The static bosonic 
Green's functions (fermionic Green's functions are never static, since the fermion frequences are 
odd) play an important role for the phase transition. For example, the expectation value of 
the scalar field is just Gi(0,0); the static correlation lengths can be extracted from two-point 
Green's functions, etc. 

Assume now that the theory is weakly coupled, and that the expectation value of the Higgs field 
in the broken phase is small enough, so that the vector boson masses are much smaller than the 
temperature. Then the description of the phase transition (i.e. static bosonic Green's functions) 
can be derived within a simpler 3d theory, which contains only bosonic fields corresponding to the 
n = sector of Fourier decomposition. In loose terms, "superheavy" (we keep the word "heavy" 
for other fields defined below) fields are integrated out. This theory is valid up to momenta 
k ^ T. To specify the dynamics of the effective theory one writes down the most general 3d 
super-renormalizable Lagrangian, containing zero modes only, and determines its parameters 
by the matching condition. This condition requires that the 2-, 3- and 4-point one-particle 



16 Whenever the comparison between 3d and 4d simulations possible, they are in agreement, indicating the 
correctness of the dimensional reduction beyond perturbation theory. Generically the errors in 4d simulations 



are considerably larger than those in 3d [ 167 1, because of rather stringent requirements on the lattice size in 4d, 
see ref. 
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irreducible Green's functions, evaluated in the effective theory and in the full 4d theory are the 
same up to some power of the coupling constant. The effective theory approximately describes 



then all static Green's fuctions of the high temperature 4d theory. As discussed in ref. ||107 
the maximum accuracy that can be reached with a super-renormalizable 3d theory is 

AC 

-^~0(g% (5.20) 

To be more explicit, take as an example the Minimal Standard Model Lagrangian. Then the 
3d effective theory is an SU(2) x U(l) bosonic theory, which contains the Higgs doublet, the 
scalar triplet (zero component of the SU{2) gauge field), and the scalar singlet (zero component 
of the U(l) field) with the action 

S = 1^x1^0^ + ^F^ + iD^^D^+ml^ + Xs^^Y + XBB^ 
+ \{D l A a Q f + \mlA%Al + \\ A {A%A%f + ^B,) 2 + ^m' 2 B B 
+h 3 &$A a A a + h' 3 &*B B - ^\ AB B &A a T a $ + XabBqAqAq }, (5.21) 

where G£ = - djAf + g 3 e ahc A\A), = - d,B u D& = (d l - ig 3 r a A?/2 + 

DiA^ = diA a Q + g 3 e abc 'A\A C Q , and $ = (0 3 + z0 4 , fa + ifa) T / ^2. The r a 's are the Pauli matrices. 

The factor 1/T multiplying the action has been scaled into the fields and the coupling constants, 

so that the fields have the dimension GeV 1//2 and the couplings gf , A3, X A , \ab have the dimension 

GeV. 

The complete one-loop calculation of the 3d coupling constants and two-loop calculation of 
the 3d masses of the effective theory can be found in ref. [ 152 , 107] ; we present here only simple 
tree relations for the coupling constants 

03 = 9 2 T, g'i = g' 2 T, A 3 = AT, h 3 = ^g 2 T , (5.22) 
h' 3 = ^g' 2 T, \ A = \ B = \ AB = 0, \ A b = gg'T 

and one-loop relations for the doublet scalar mass 

m 2 M = -\m 2 H + tQa 3 + ^g 2 3 + ^g' 3 2 + \qy) (5-23) 



J2 



m D = [- + 



1 , 5u f\ ,2^2 



6 9 



)g' 2 T\ (5.24) 



and triplet and singlet masses 



Here m H is the zero temperature Higgs mass, gy is the Yukawa coupling constant corresponding 
to the t-quark, np = 3 is the number of fermion generations. 
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Clearly, is is much easier to study this 3d model than the initial 4d theory, since it contains 
much less degrees of freedom. Another important fact is that the only remnant of fermions is 
their contribution to the 3d masses and coupling constants. 

An inspection of relations ( 5.23[ ) and ( |5.25| ) shows that further simplification of the effective 
theory is possible. Indeed, at the tree level at m\ < the symmetry is broken, and at m\ > 
it is restored. In other words, the phase transition takes place near m\ = 0. At this point the 
hierarchy m\ <C m 2 D and m\ <C m% holds, allowing for the construction of the effective theory 
containing the Higgs doublet and SU{2) x £7(1) gauge fields only. The "heavy" scale ~ gT is 
integrated out, and the super- renormalizable Lagrangian of the effective theory is merely 

L = \ G % G % + \ F a F a + (A$ ) t (A$) + m^ f $ + A 3 ($ f $ ) 2 . (5.26) 
In the one-loop approximation 

mg(/i) = m\ — — (sh^mD + h'^m'jj] (5-27) 

and A3 = A3 at the tree level. Higher order corrections to these relations can be found in ref. 
152], p.07|| . In fact, the possibility of integrating out the heavy scale has a general character and 



can be applied to any gauge theory. The effective description in terms of the super-renormalizable 
Lagrangian for the "light" modes only provides an accuracy 

AC 

— ~0(g s ). (5.28) 

for Green's functions for "light" (static gauge bosons and scalars) fields. This theory is valid for 
k <C gT, but k may be as large as g 2 T. The effective potential for the scalar field in 3d coincides 
with the hard loop resummed potential in 4d at high temperatures ||149| , |151] , |152|| . 



Since all the four parameters of the 3d theory are dimensionful, the theory is uniquely fixed 
by three dimensionless ratios: 

y=% (5.29) 

93 93 \93) 

and the overall scale g 2 , on which the physics does not depend. The temperature variation 
changes the parameter y, while x and z depend on T only through logarithmic corrections. So, 
we arrive at the important conclusion that the dynamics of the electroweak phase transition 
depends on one dimensionless number only (x) since z is related to the known weak mixing 
angle, z ~ 0.3, and y is fixed by the requirement that the temperature is equal to the critical 
one. The dependence of the parameter x on the physical Higgs mass in MSM was found in ref. 
T07| , see fig. |. 



In fact, the 3d theory defined by ( |5.26| ) plays the role of the universal theory describing 
the phase transition not only in the Minimal Standard Model, but in many of its extensions, 
including sypersymmetry, two Higgs doublets, etc. [ [I07| ]. Let us take as an example the two 
Higgs doublet model. The integration over the superheavy modes gives a 3d SU{2) x U{1) 
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theory with an extra Higgs doublet in addition to the theory considered above. Construct now 
the one-loop scalar mass matrix for the doublets and find the temperatures at which one of its 
eigenvalues is zero. Take the higher temperature; this is the temperature near which the phase 
transition takes place. Determine the mass of the other scalar at this temperature. Generally, 
it is of the order of gT, and therefore, is heavy. Integrate this heavy scalar out together with 
the heavy triplet and singlet - the result is the simple SU(2) x U(l) model. In the case when 
both scalars are light near the critical temperature, a more complicated model, containing two 
scalar doublets, should be studied. It is clear, however, that this case requires fine tuning. The 
consideration of the phase transitions in the two Higgs doublet model at the one-loop level can 
be found in refs. JT8|, |8|, H HM [L87l |8|, |18|, |190| . 

The same strategy is applicable to the Minimal Supersymmetric Standard Model (MSSM). If 
there is no breaking of colour and charge at high temperature (breaking is possible, in principle, 
since the theory contains squarks), then all degrees of freedom, excluding those belonging to the 
two Higgs doublet model, can be integrated out. We then return back to the case considered 
previously. The conclusion in this case is similar to the previous one, namely that the phase 
transition in the MSSM can be described by a 3d SU(2) x U{1) gauge-Higgs model, at least in 
a considerable part of the parameter space. A one-loop analysis of this theory was carried out 



in refs. [191, 192, 193, 191 



Of course, effective super-renormalizable 3d theories do not provide an exact description of 
the phase transitions. Parametrically, ~ 0(g 4 ) or ~ 0(g 3 ) depending on the level of dimen- 
sional reduction. The numerical estimate of uncertainty is model dependent. For the standard 
electroweak theory, corresponding estimates were given in ref. ||107|| , with the result that the 
effective description provides an accuracy ~ 1 % for Higgs masses from 30 GeV up to few hundred 
GeV. For smaller Higgs masses the phase transition is too strong and the vector boson masses 
are of the order of temperature, i.e. the assumption of the scale hierarchy does not hold; if the 
Higgs mass is close to the unitarity bound, the perturbation theory, used in the construction of 
the effective Lagrangian, breaks down. 

To summarize this discussion, in order to study phase transitions in weakly coupled gauge 
theories, one may construct simpler 3d theories, parameters of which can be perturbatively 
calculated. This calculation is free from any infrared divergences, and its ultraviolet divergences 
are removed by the usual counterterms of the zero-temperature perturbation theory. A unique 
3d theory plays a role of universal theory for the description of the phase transition in many 4d 
models. In particular, the strength of the electroweak phase transition depends on the unique 
number, x = Xz/dl- 



5.4 Phase structure of 3d theory 

The construction of the effective 3d theory is free from infrared divergences. However, the 
perturbative calculations in 3d are infrared divergent in the symmetric phase. This is easy to 
see by simple power counting. For example, the iV-th loop contribution to the vacuum energy in 
3d (with dimension GeV 3 ) must contain the factor (gl)^ -1 . Therefore, starting from N = 5 some 
dimensionful parameter (say, the typical momentum scale k) should appear in the denominator 
of each term of the expansion, i.e. the expansion parameter is g\jh. In the symmetric phase 
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the infrared cutoff is absent, and perturbation theory breaks down. The first (logarithmically) 
divergent contribution in the vacuum energy appears in four loops and is of the order of (gf) 3 , 
giving a contribution of the order of g 6 T A to the free energy. 

This means that the 3d theory should be treated in a non-perturbative way. At present, 
non-perturbative studies of the complete 3d 577(2) x £7(1) model are absent, and most results 
were derived for an SU (2) gauge-Higgs system. Let us consider first the general structure of the 
phase diagram of this theory. In the SU(2) gauge-Higgs model there are no known local gauge 
invariant order parameters, which would acquire a non-zero value only in one of the two phases 
of the model. In other words, there is no symmetry breaking or restoration in this model, and 
the gauge symmetry is always intact (see, for instance refs. [|195| , |196| ). An ordinary jargon, 
saying that the vector bosons are massless in the symmetric phase and massive in the broken 
one is incorrect, since vector degrees of freedom are massive at any choice of parameters. A 
good statistical analogue of this system is an ordinary liquid, which may undergo the first order 
liquid-vapour phase transition without restoration or breaking of any symmetry. The phase 
diagram of the latter system usually has a critical point, where the line of the first order phase 
transition ends. At this point the phase transition is of the second order, and there is a massless 
scalar excitation in the physical spectrum. 

The absence of the true order parameter in the gauge-Higgs system suggests two possible 
phase diagrams in the (x, y) plane, shown schematically in fig. ||. In the first case the line 
of the first order phase transitions has an end point at finite x, while in the second case the 
end point is at x = oo. By dotted lines we show the upper and lower spinodial decomposition 
temperatures. In the first case there exists a critical scalar self-coupling above which there is 
no phase transition at all, while in the second case a first order phase transition occurs for all 
physical parameters. It is important for cosmological applications to clarify the type of the phase 
diagram. For example, if there exists an end point of the line of the first order phase transition, 
and the physical value of x is to the right of this end point, then there is a smooth cross-over 
transition without any strong deviations from thermal equilibrium. In the latter case there are 
no observational consequences from the electroweak epoch. 



The following approaches were applied to SU (2) gauge-Higgs system: e expansion ||197|| 



exact renormalization group approach [|198| , |199| , |200| , |201|| , Schwinger-Dyson equation ||202 



perturbation theory | |152| |, and lattice Monte Carlo simulations | ]169| , |203| , |204| | . The first two 
approaches favour the second type of the phase diagram, while the third one argues in favour 
of the end point of the critical line. In general, perturbation theory describes well different 
characteristics of the phase transition at small scalar self-coupling constants, while e-expansion 
is more suitable in the regime when X^/gl ^ 1 |205| , |206| , |207| |. The lattice MC simulations 
cannot resolve the order of the phase transition at large scalar self-coupling constant. 

We shall review here the results from the lattice 3d MC simulations only [[208| , |209| , |168| , p.69| , 



203], |210|| . The comprehensive presentation of all existing approaches would make the discussion 
too lengthy; our emphasis on MC simulations is mostly because this is the only approach which 
is based on the first principles and does not require any extra assumptions. For example, e 
expansion relies on the hope that e = 1 may appear to be a small expansion parameter. Exact 
renormalization group approach requires some truncation of the equations, and Schwinger-Dyson 
equations are constructed with the use of the perturbation theory which may appear to be not 
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applicable. The lattice MC simulations produce a set of "true" characteristics of the phase 
transition (provided the quality of the data is such that extrapolation to the continuum limit 
is possible), which then can be used in cosmological applications. Of course, MC simulations 
do not provide an analytical understanding of the dynamics of the transition, but they may 
be considered as an "experimental" basis for construction or testing of different theories of 
the phase transitions. The lattice MC simulations fail to describe very weakly first order phase 
transitions. However, only sufficiently strong first order electroweak phase transition is of interest 
for cosmology, i.e. the lattice simulations can cover the entire interesting range of parameters. 

5.5 The lattice formulation 

The 3d gauge theory defined by the action 

S = J (fx^Gtptj + (A$) f (A$) + rn\&§ + A 3 ($ f $) 2 } 



(5.30) 



has a number of remarkable properties. It contains only dimensionful coupling constants, and, 
therefore, it is super-renormalizable. The only ultraviolet infinities are those in the scalar mass 
renormalization, while the (3 functions for g\ and A3 are equal to zero. The exact (3 function for 
the mass parameter is known ||152|| . 
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h m = tttSs + 9A 3 <7 3 2 - 12A*. (5.32) 
lb 

The lattice version of the continuum theory ( |5.30p is defined by the following action: 

x i<j - ' 

- &EE^ $t WW( x + l )+ (5-33) 

x i 1 

+ E^r$t(x)$(x) +^E[^Tr$ t (x)$(x) - l] 2 . 

X Z X Z 

Here Ui(x) and are the link and plaquette variables, and $ is the scalar field. The action 
depends on the three dimensionless parameters Pg:Ph,Pr (recall that the continuum dynamics 
is completely described by two numbers, x and y) . The fact that the theory under consideration 
is super-renormalizable enables one to find an exact (in the continuum limit) matching between 
the lattice variables and the physical parameters. The corresponding calculation has been carried 



out in refs. ||168| , |211|| ; it provides the relation between the lattice renormalization scheme and 
the MS scheme in continuum. The continuum limit is (3q — > 00, (3 H — > ~, /3r — > in such a way 



35 



that the physical parameters x and y, defined as 



1, PrP, 



x = ^-X 3 a(3 G = (5.34) 
4 Ph 



^(±-3-^1 + ^(1 + 4*) + 
8 \p H (3 G ) + 327T ^ i + 4 ^ + 

+ 9a; - 12x 2 ) fin ^ + A + V + V x 



+ 
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(5.35) 



stay constant. Equation ( |5.35| ) depends on several constants arising from lattice perturbation 
theory; £ = 3.17591, ( = 0.09, r\ = 5.0 and f) = 5.2 were computed in refs. ||168| , [211|| . 



If we denote by a the lattice spacing which has dimension of length, then the overall continuum 
scale g\ is related to it as 

4 

g\a = — , (5.36) 

Pg 

i.e. the continuum limit corresponds to a — > 0. The relation of the lattice variables to the 
continuum fields in this limit is 

$ = VR L , R 2 L = j-$4> = ^Tr Ui = exp (^iag 3 T a A^ , (5.37) 

where V is a unitary SU(2) matrix, Rl is the radial mode of the Higgs field, r a are the Pauli 
matrices. 

Equations (|5.34j) and ( ^.35| ) serve as the basis for relating the results of the lattice simulations 
to the physical reality. We are not going to discuss here technical details of the lattice simulations 



169|| . Instead, we review the general strategy for the lattice study of the phase transition. The 
reader not familiar with MC computations may consider the computer as a device for calculating 
the Euclidean functional integral 

{0{A 0)) = JDAD^eMS) ' (5 ' 38) 

where 0(A, (j>) is some gauge-invariant functional of the gauge and scalar fields, and S is the 
lattice action. The integral should be computed for different volumes of the system and different 
lattice spacings, and extrapolation of the results to the infinite volume and zero a must be taken 
at the end. By choosing the specific 5-function form of this functional, 

O = 5{c-G(A,<f))) (5.39) 

one can construct probability distributions for various order parameters G(A,<p), essential for 
the study of the phase transition. 

The critical temperature. Let us fix /3q (i.e. the ultraviolet cutoff) and vary parameters 
(3h and (3r in such a way that the physical variable x stays constant. In 4d language this 
corresponds to varying the temperature. The existence of the first order phase transition means 
that at some critical value of /3# there are two degenerate ground states with different properties. 
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In particular, expectation values of the gauge invariant observables are different in distinct states. 
For example, the average value of R 2 is expected to be smaller in the symmetric phase than in 
the Higgs phase. This means that the probability distribution at the critical point has non- 
trivial form and contains two peaks rather than one. The typical picture of the distribution 
evolution is shown in fig. |5|. At large fin (low temperatures) there is a unique value of the 
order parameter R 2 , corresponding to the Higgs phase, while at some (3h a double peak stucture 
develops, signalling about the first order phase transition. The critical value of (3h can be found 
from the requirement that the areas under the two peaks are equal to each other. Of course, 
it must be checked that the double peak structure is not a lattice artefact and passes different 
tests singling out the first order behaviour |169| . After the critical value of (3h is found, it may 



be converted to the critical temperature of the underlying 4d theory. 

The latent heat. The latent heat L — the energy released in the transition — can be 
calculated from 

L dAp T d . , „ T d . „ ,„ 

- = — - = A log Z = AP , 5.40 

T dT VdT & VdT K J 

where the derivatives are evaluated at the critical temperature, Ap is the difference of the 

pressures in the symmetric and broken phases, and AP is the difference in the probabilities of 

the phases in volume V. In eq. (|5.40|) , T is the physical (4d) temperature. The quantity AP is 

directly proportional to the difference of the areas under the two peaks in the order parameter 

distributions near T c ; it may be computed by performing simulations at (3h close to the critical 

one. 

The jump of the order parameter. The order parameter usually discussed in the study 
of the effective potential is the vacuum expectation value v of the Higgs field. This quantity, 
however, is not gauge invariant. The gauge-invariant substitute of it is the scalar condensate 
(0V)- This is a composite operator, whose expectation value contains linear divergence at 
the one-loop level and logarithmic divergence at the two-loop level. These divergences can be 
removed, e.g. by the MS prescription; the resulting condensate is then dependent on the scale 
parameter /x. Because of the fact that the 3d theory is super-renormalizable, an exact relation 
between the lattice quantity (R 2 ) and {$4>) can be found ||168|| , 



M ~ \ Mb Ur>) ^g ^1 + C + W - *)+ O (±- ). (5.41) 



gl s yuyn V 1 ttPhJ (4tt) 2 V ta 2n s 4 J VA 



>G . 



Numerically £ + ^S 2 — 5 = 0.67. Thus, the extrapolation of lattice measurements of the quantity 
(R 2 ) to the limit (3q ~^ 00 allows one to determine an "exact" value of the scalar condensate. 
The comparison of the lattice results with known two-loop perturbative expansion allows the 
extraction of the magnitude of the higher order perturbative terms. In this way a three-loop 



correction to the effective potential has been numerically determined in ref . ||169|| . For estimates 



the relation \%a = (ft(f)(T))/T can be used. The effective potential for the gauge-invariant 



condensate was constructed in ref. 212 



The jump of the quantity ((jrcfi) at the phase transition (difference between condensates in 
the broken and symmetric phases at the critical temperature) is a finite and scale independent 
quantity, which can be directly measured on the lattice as the distance between the positions of 
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the two peaks in the R 2 distribution, 



A((0V)) = - g 2 p G [3 H A((R 2 )). (5.42) 
It can be shown that A ((0^0)) is directly related to the latent heat of the transition [|168|, |169| . 



213 



The interface tension. The interface tension is one of the most important quantities which 
characterize the strength of the phase transition. It can be measured by constructing probability 
distributions of some local order parameter. At the critical temperature, a system in finite volume 
predominantly resides in either the Higgs or the symmetric phase, but it can also exist in a mixed 
state consisting of domains of the two states. The probability of the mixed state is suppressed 
by the extra free energy associated with the interfaces between the phases. This causes the 
typical two-peak structure of the probability distribution of the order parameter at the critical 
temperature (see fig. |G|): the midpoint between the peaks corresponds to a state which consists of 
equal volumes of the symmetric and broken phases. Because of the associated extra free energy, 
the area of the interfaces tends to minimize. Assuming lattice with periodic boundary conditions 
and geometry L 2 x L z , where L x < L z , the minimum area is 2 x A = 2(L x a) 2 — the factor 2 
appearing because there are two separate interfaces. The interface tension a can be extracted 
from the limit ^ 

1 v^oo 2A R min 

where P max is the maximum of the probability distribution in the peak and P m i n is the minimum of 
distribution between the peaks in fig. At sufficiently large size of the system in the z direction, 
the probability distribution has a characteristic plateau, corresponding to the translational zero 
mode of the domain wall. Other methods for determination of the surface tension are discussed 
166||, the estimate of the higher order perturbative corrections is contained in ref. |21 \\ . 
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5.6 Some lattice results 

The 3d lattice simulations have been made for four values of the continuum parameter x, namely, 



x = 0.01830, 0.06444, 0.08970 and x = 0.1188 p9l ETA. These values correspond to different 



4d physical parameters in different models. For definiteness, let us take the SU(2) sector of the 
standard model with top quark mass m t = 175 GeV. Then the first value of x cannot be realized 
with any value of the Higgs mass (see fig. §), while the others correspond to inn ~ 51.2, 68.0 
GeV and ^ 81 GeV, respectively. The latter numbers come from the one-loop relations between 
the MS and physical parameters of the standard model and the one-loop dimensional reduction 



107|1 . Because of the large value of the top Yukawa coupling constant the one- loop corrections 
are quite substantial for small Higgs masses, for example Smn/friH ~ 15% for m# »s 50 GeV. A 
naive estimate of the two-loop corrections (taken as square of the one-loop contribution) indicates 
that the accuracy of the one-loop computation of the physical masses at fixed 3d parameter x 
is ~ 2% for x = 0.06444 and better for larger x. 

For the "large" value of the Higgs mass (81 GeV) it was not possible to resolve the order of 
the phase transition on the lattices up to 48 3 , and the data are compatible with the smooth 
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cross-over, second order phase transition or very weakly first order phase transition. For smaller 
Higgs masses the transition is of the first order. 

The most complete study has been done for x = 0.06444. This value may be quite realistic 
for extended versions of the electroweak theory, but is excluded experimentally in the MSM (it 
corresponds to m# ~ 51.2 GeV). Below we present the specific numbers for MSM with m t = 175 
GeV, derived with the use of one-loop relations. The critical temperature of the phase transition 
is T c = 89.79 GeV, and the vev-to-temperature ratio at T c is v/T c = 0.64. The domain wall 
separating the broken and symmetric phases has the surface tension o ~ 0.002T C 3 , and the latent 
heat is L/T^ = 0.12. At T > T + = 89.93 GeV only the symmetric phase is stable, and at 
T < 71 = 89.36 only the broken phase is stable, while at temperatures between T_ and T + 
both phases can exist simultaneously. The scalar correlation lengths in the symmetric and Higgs 
phases are ~ 6/T c and ~ 8/T c , respectively. The statistical errors of the lattice numbers are 
0.015% for the critical temperature, 1% for the expectation value of the Higgs field and the 
latent heat, and ~ 20% for the interface tension. The two- loop corrections may introduce extra 
uncertainties in these numbers of the order of 2%. 

The bubble nucleation temperature T bubb i e lies somewhere between T c and T_ and may be 
estimated with the use of the surface tension and latent heat found on the lattice. Inserting 
the lattice numbers to the relations ( |5.5| ), ( |5.8|) , ( |5.7| ) gives an estimate ~ 0.0004, i.e. the 
bubble nucleation temperature is very close to the critical one. The smallness of AT/T C is due 
to fact that the ratio cr 3 /L 2 T c ~ 2 • 10~ 6 is so small [gT§. Since AT/(T C - H) ~ 0.1 is also 
small, one is in the thin wall regime; indeed, the size of the bubbles when they nucleate is at 
least R c ~ 110/T C , which is much larger than the scalar correlation length in either the broken 
or the symmetric phase at T bubble . Because T bubble is very close to the critical temperature, the 
expectation value of the Higgs field at T bubb i e is almost the same as at T c . In other words, the 
transition is very weakly first order. 

The lattice results can be compared with the perturbative calculations of the two-loop effective 
potential. The predictions of the critical temperature, latent heat, and the jump of the order 
parameter appear to be quite reasonable and are within a few per cent for these quantities. 
However, perturbation theory fails to describe the surface tension (and, therefore, the bubble 
nucleation rate) at least for x = 0.06444: the perturbative value is about 3 times larger than the 
lattice one. 



5.7 Dynamics of the phase transition 

In the case of the "vacuum phase transitions" - false vacuum decay at zero temperatures - 
the energy stored in the metastable vacuum transforms into kinetic energy of the domain walls 
126| , |127| |. As a result, the velocity of a domain wall increases and approaches the speed of light. 



Then in the collisions of domain walls, their kinetic energy is released and transforms into heat, 
and the system is reheated up to a certain temperature. 

At non-zero temperatures nucleated bubbles expand in the medium, and the latent heat of 
the transition may be released by many different mechanisms. In the idealistic case of the very 
slow expansion of the Universe, the temperature stays constant during the phase transition, fig. 
[7| (see, e.g. ref. ||138|| ). Every point on the plateau corresponds to a mixed state containing 
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domains of broken and symmetric phase; the left point of the plateau corresponds to the pure 
symmetric phase, and the right point to the pure broken phase. The release of the latent heat 
of the transition is accommodated by the Universe expansion. 

In fact, the Universe expands not so slowly at the electroweak scale and is supercooled in 
the symmetric phase. One can distinguish several different epochs in the phase transition. The 
first one is the bubble nucleation, the second one is the bubble growth, the third one is the 
bubble percolation - the period when different bubbles collide. At this stage the Universe may 
be reheated up to the critical temperature (fig. 0). If it happens, then the later evolution may 
be close to an ideal case described above. 

We have already discussed the bubble nucleation rate and determination of the bubble nu- 
cleation temperature Let us now consider the bubble expansion from macroscopic point of 
view in more detail |17|, |138|, |T|, |5|, gT5j pOLlSkj. 

Initially, the bubble of a new phase is a microscopic object with a size of several correlation 
lengths. At the bubble nucleation temperature the energy density in the symmetric phase is 
larger than that in the broken phase, 
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(5.44) 



Here N e ff = Nb + |n/ is the effective number of the massless degrees of freedom and L is the 
latent heat. For the Minimal Standard Model N e ff = 106|. The general hydrodynamical con- 
sideration of the bubble evolution leads to two possible types of bubbles, known as deflagration 
and detonation bubbles 



2T6|, gr?|. 



Consider an isolated macroscopic bubble. For a large enough bubble, its curvature may be 
neglected and the interface may be taken as a planar domain wall Let us proceed in the rest 
frame of the domain wall and denote by v i the velocity of the medium falling on the domain wall 
(symmetric phase) and by v 2 the velocity of the medium inside the bubble. Then, if v\ < c s , 
where c s ~ ^= is the velocity of sound in the medium, then the phase transition proceeds 
through deflagration. In this case the medium is accelerated when it passes through the domain 
wall, t>2 > v i- If; on the contrary, V\ > c s , then v% > v 2 , and we have detonation bubbles. 
The realization of one of the two mechanisms of the bubble walls propagation depends on the 
relationship between the latent heat of the transition, surface energy density, and the rate of 
entropy generation on the phase boundary ||217|| . The latter is to be found from the microscopic 
analysis of the interaction of particles with the bubble walls ||223. 



If we now choose the reference frame to be the rest frame of the plasma before the bubble has 
nucleated, then the medium should be at rest in the centre of the bubble and far from it. Then, 
for the deflagration bubbles the medium in front of the bubble wall is accelerated by the motion of 
the wall, and there is a shock wave in the symmetric phase moving with the velocity v s h oc k > Vb, 
and Vb = v 2 is the bubble wall velocity. The velocity of the symmetric phase plasma between 
the shock front and the bubble wall is given by v sym = {v 2 — fi)/(l — viv 2 ). The temperature 
of the medium between the shock front and the bubble wall Tg^ock an d the temperature inside 
the bubble T r are different from the bubble nucleation temperature (temperature outside the 



17 As was shown in 



the domain walls are stable against small perturbations. 
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shock front) They depend in general on the distance from the bubble centre. Usually the 
inequality T& < T r < T s h oc k holds true ||217||, but T r < T at high rates of entropy generation at 



the domain walls ||215|| . The latent heat of the transition transforms into the kinetic energy of 
the plasma in the symmetric phase, and to heating of the plasma inside the shock wave front. 

The detonation bubbles have a different structure. The velocity of the domain wall is larger 
than the speed of sound, the symmetric phase plasma is at rest right in front of the domain wall 
and has temperature T&. The plasma just behind the wall (in the Higgs phase) is accelerated 
by it and has velocity v br = (v\ — i>2)/(l ~ Wa)- Finally, the plasma is stopped at some surface 
inside the bubble by a rarefaction wave. As in the previous case, the temperatures inside the 
front of the rarefaction wave T r and behind the bubble wall T& r are different from T&. 

The recent analysis carried out in ref. ||22£ , 227 | suggests that for a sufficiently wide range 
of the parameter space of the standard model, m H < 90 GeV, an isolated bubble expands as a 
weak deflagration ("weak" means that the velocity of the bubble wall is subsonic). The velocity 
Vb was found to be in the range 0.38 < v & < 0.45, while the velocity of the shock front is close to 
the speed of sound in the medium, v s h oc k — Cs — 0.58. When the shocks coming from different 
bubbles begin to collide, the single bubble approximation breaks down. At this time, roughly, 
a (vb/c s ) 3 ~ 0.3 part of the Universe volume is in the broken phase. The subsequent evolution 
of bubbles depends on the temperature in the symmetric phase, which will emerge as a result 
of shock waves interaction. This temperature can be estimated as follows [ 215 , [229|| . Suppose 
that all latent heat of the transition is immediately released in form of heat in the broken phase. 
Then the reheating temperature can be found from energy conservation condition, 

2 2 

h Nef,T > = h NeffT " ~ L ' (5 - 45) 



so that 

T r -T h 15 L 
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(5.46) 



Now, if Tc ~ Tb 3> Tr y Tb , then the reheating process may be neglected. However, in the opposite 
case the heat release is important and the expansion of the bubbles of a new phase should slow 
down. 

Simple estimates can be done in a thin wall approximation. The bubble nucleation temperature 
is given by 

T c ~T b / 16trx 3 

-fr = izm^' (5 - 47) 

where Sq ~ 160 is the action for the critical bubble and a is the surface tension. For example, 
for the Minimal Standard Model with m# = 51.2 GeV and m t = 175 GeV we get T c = 88.93 
GeV, L/T A = 0.124, a/T 3 = 0.0023 pm and 



Tl—It = 9 • 10- 4 > ^— = 3 ■ lO" 4 . (5.4? 



18 Of course, no such Higgs boson exist in minimal standard model because of the experimental constraint. 
This value is taken because this is the highest MSM Higgs mass for which the magnitude of the surface tension 
is known reliably from the lattice simulations, see section 5.5. Being unrealistic for MSM, this example is 
phenomenologically acceptable for the extensions of the standard model. 
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So, for this choice of parameters the Universe should reheat to the critical temperature, the 
speed of the domain walls is greatly reduced, the broken phase bubbles expand slowly due to the 
Universe expansion. At the final stage of the phase transition the remnants of the symmetric 
phase shrink, again due to the Universe expansion. This picture is also true for higher values of 
the MSM Higgs masses, provided the phase transition is still of the first order. 

An estimate of the bubble wall velocity at a late stage of the phase transition (under the 
assumption that the Universe is reheated up to the critical temperature) can be found from 
a simple thermodynamical consideration |230|1 (see also ref. ||229| ). Suppose that the average 



bubble size is Rb, and that the part of the volume of the space occupied by the broken phase is 
V. Then the requirement that the Universe expands adiabatically is given by 

^VL = sH, (5.49) 
Rb 

where s is the entropy density in the symmetric phase, s = ^-N e ffT A , and H is the Hubble 
constant. The average bubble size can be found from the following consideration. The (non- 
normalized) distribution of the bubbles in sizes at moment t is given by 

P{R,t)dR~1l{h), (5.50) 

where 1Z(ti) is the bubble nucleation probability at time t\, and ti is related to R and t by the 
obvious condition that U = t — — . Then 



Vb 



J R P(R,t) dR , . 

{R){t) = J P(R,t) dR ■ (5 - 51) 

In the thin wall approximation 

with A = §|^fr, and t c is the time corresponding to the temperature T c . The average bubble 
radius at time t is 

(R)(t) ~ -v K -—^- (5.53) 



2 At\ 

and at the percolation temperature it is about (R)(t) ~ l^c-y^f • For example, for the numerical 
values of the parameters given above one finds A ~ 5.3 • 1CT 5 , HRb ~ 1CT 6 and 

iU^^^-KT 3 (5.54) 

for V ~ 0.3. This value gets larger if the scalar self-coupling is decreased. For example, for 
stronger phase transition the reheating up to the critical temperature may occur, and the velocity 
of the domain wall from the estimate of ref. [ [229|| may be higher by a factor ~ 10. 



Of course, the estimates given above are rather rough since they are based on the assumption 
of the instantaneous latent heat release. Nevertheless, they show that it is quite plausible that 
the slow stage of the phase transition takes place for the interesting range of the parameters of 
the underlying theory ||229| . 
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6 Survival of primordial baryon asymmetry 



The anomalous electroweak processes are rapid at sufficiently high temperatures. Their rate 
T sp h exceeds the rate of the Universe expansion |^ in the standard Big Bang scenario in the 
following interval of temperatures: 

lOOGeV ~ T* < T < T** ~ a^M Pl ~ 10 12 GeV, (6.1) 

where the lower temperature T* is to be found from the condition of decoupling of the sphaleron 
processes in the broken phase of the EW theory [1231 , 232], 



y * 45. (6.2) 

Clearly, the equilibrium character of 5-violating reactions has an important impact on the 
survival of the primordial baryon asymmetry. Several different cases can be distinguished, de- 
pending on initial conditions and on the rate of B and L non-conservation due to processes other 
than those associated with sphalerons. 

(i) Suppose that the Universe is asymmetric with respect to the anomaly free fermionic charges 
A,- = La — —B of the standard model at T > T**, and assume that at T < T** there is no B or L 
violating interactions besides the electroweak anomalous processes. The origin of the primordial 
asymmetry is not essential here. Then anomalous reactions convert the initial asymmetry to the 
baryonic one at T = T*. For the minimal standard model the relationship is given by |94], [H| 
(see also @, |7|), 



A - 8n f + 4 A K 4 V m ^ A (6 3) 

A ° " 22n f + 13 Ab - L " ^13^ h l^f^ (6 ' 3) 

where mf is the lepton mass of a given generation, K « 1. The first term in the right hand side 
of this equation tells that (B — L) asymmetry is reprocessed into the baryon asymmetry, while 
(B + L) tends to be washed out; the second term is the correction coming from slightly different 
behavior of quarks with different masses in the plasma. If the initial value of (B — L) is non-zero 
(coming, say, from the GUT physics) then the baryon asymmetry, up to a possible contribution 
from the EW phase transition (see below), has a primordial character. If, on the contrary, the 
initial B — L asymmetry is absent, we can rely only on the second term in ( |6.3j ). For three lepton 
generations one gets a suppression A — 3 x 1CT 6 A 3 . So, to have a non-negligible effect, the 
initial asymmetry A3 must be very large, or the standard theory should be extended by adding 
heavy leptons. 

(ii) Suppose now that there are some reactions, which do not conserve all Aj, and which are in 
thermal equilibrium for some period between T* and T**. At this intermediate epoch B and L are 
non-conserved separately, and according to the third Sakharov condition all baryonic and leptonic 
asymmetries are washed out. Hence, the existence of these reactions is fatal for the primordial 
baryon asymmetry. If the baryon asymmetry is not produced at a later time, the requirement of 
the absence of these reactions may appear to be a powerful tool for constraining the properties of 



new particle interactions 233, 96, 234, 235 . However, some time ago it was realized that most 
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of these constraints are drastically weakened due to the smallness of some Yukawa coupling 
constants in the standard model or its supersymmetric extensions ||236| , [237| , |238| , [239| , |240| , |241| . 

Let us discuss the main idea of these estimates on the example of lepton number violating 
interactions, leading to the Majorana neutrino masses [ |240| . We take for simplicity the 
Minimal Standard Model and add to it lepton number violating interactions. The SU(2) x U(l) 
symmetric low energy Lagrangian with AL = 2 has the form 



1 



(6.4) 



where Li and L°- are lepton doublet and its charge conjugate, respectively, <fi is the scalar doublet, 
v is the vacuum expectation value of the Higgs field v = 246 GeV. The rate of L non-conserving 
reactions Lip — > L c <fi* at high temperatures has been found in ref. ||240 



r 



9 T 3 n 

7T b V 4 



(6.5) 



where m v is an average Majorana neutrino mass, 



m, 



I I 2 1 I |2 1 I |2 

\m ee \ + |ra eAt | + \m eT \ . 



(6.6) 



These reactions were initially required to be out of thermal equilibrium at T < T**; this 
leads to a very stringent constraint 



m v < 



a/M T** 



l(T 2 eV. 



(6.7) 



An implicit assumption in the derivation above is that the set of conserved numbers Aj is 
a complete one below T = T**. In fact, this is not true due to the smallness of the right- 
handed electron Yukawa coupling constant. In the limit when this constant is zero, the right- 
handed electron number is conserved, and the asymmetry in it propagates to the asymmetry in 
baryon number. The rate of reactions not conserving the right-handed electron number (say, 
eiH — > ep>W) is of the order of 

Tr ~ ccwflT (6.8) 

where f e is the electron Yukawa coupling constant. These reactions are out of equilibrium at 
T > T R ~ 3 TeV | |240| | . It is this temperature which should be used in eq. (p-7\) instead of T** . 
Thus, we arrive at much weaker constraint |[239| , |240| , |2 4 2 



m v < 8KeV. 



(6.9) 



which must be satisfied in any case because of the known laboratory limits and other cosmological 
considerations. 

The same type of considerations apply to other possible interactions breaking lepton and 
baryon numbers. The general conclusion is that the initial charge asymmetry can survive during 



44 



the epoch at which anomalous reactions are at thermal equilibrium. Moreover, initial asymme- 
tries in fermionic quantum numbers, different from the baryon number, are usually transferred 
to baryon asymmetry towards the end of the equilibrium sphaleron period. 

We barely know the history of the Universe at very high temperatures (say, at T ^> 1 TeV). 
It may well be that the Universe was symmetric with respect to all fermion charges at T > 10 12 
GeV. This assumption, being a bit arbitrary, may be in fact a natural consequence of inflation, 
which exponentially dilutes the densities of all global quantum numbers (e.g. baryonic or lep- 
tonic). If true, baryon asymmetry should be produced at relatively late stages of the Universe 
expansion. As pointed out in section 3, this may happen either at intermediate temperatures (1 
TeV < T < 10 12 GeV) or at the electroweak temperature (T ~ (a few) x 100 GeV). 

Our main topic is the discussion of baryogenesis in the case where the only relevant source of 
B and L non-conservation is the electroweak anomaly. In a sense, this is the most conservative 
possibility, since it relies only on physics we trust experimentally. We further constrain ourselves 
and consider the Minimal Standard Model or its natural extensions, such as the two Higgs 
doublet model or supersymmetry. In these models the only known possibility to generate the 
observed baryon asymmetry is that associated with the electroweak phase transition. Further 
extensions of the standard theory, containing topological defects (such as strings) can also lead 
to baryogenesis via anomalous reactions. The discussion of this interesting possibility can be 
found in refs. pi |44|, B HI H3- 



7 Electroweak baryogenesis 
7.1 Strength of the phase transition 

There is a general condition that should be satisfied in any particle physics model used for the 
generation of the baryon asymmetry at the electroweak scale [ [231 , [232| , |§5| . Namely, the baryon 



asymmetry created by some mechanism must not be erased by the anomalous reactions. In 
other words, the sphaleron processes should be out of thermal equilibrium immediately after 
the electroweak phase transition (in the broken phase), i.e. inequality ( |6.2| ) must be satisfied. 
The temperature T* can be as large as the critical temperature T c , if the Universe is reheated 
up to it, or as small as the bubble nucleation temperature T b . The requirement ( |6.2| ) places a 
strong constraint on the strength of the phase transition, and, therefore, on the parameters of the 
electroweak theory. The recent discussion of this bound, incorporating the results of numerical 



simulations of the electroweak phase transition is contained in ref. ||169|| . We sketch here the 
main points of the analysis. 

As discussed in Section 5, the 3d SU(2) x U(l) gauge-Higgs theory plays a role of universal 
theory of the electroweak phase transition in the Minimal Standard Model and a number of 
its extensions. The 3d effective theory is characterized by a unique parameter, x = Xs/g^, 
completely defining its dynamics; in particular, the effective sphaleron mass is a function of it. 



In the one-loop approximation ||109| , |110| 



E sph (T) D (X 3 \ 
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where is the scalar field expectation value determined from the one-loop effective potential. 
The two-loop corrections to the sphaleron mass are unknown; parametrically 



E sph \nm T J 

where ttit is defined by eq. ( |4.40| ) . The perturbative and numerical analysis of various quantities 
in the broken phase (such as free energy, correlation lengths, vacuum expectation value of the 
Higgs field) suggests that the "true" expansion parameter is Kgl/{nm T ) with k ~ 1. So, it is 



natural to assume that \A\ ~ 1. Then, from eqs. (^2) and (7\1) one obtains v/T > 1.2 or 
v/T > 1.5, depending on the sign of A. To establish a conservative upper bound, a bubble 
nucleation temperature T&, which is somewhat smaller than the critical temperature, should 
be taken. If the perturbative description of the bubble nucleation based on two-loop effective 
potential is valid, then at T& the ratio v/T is about 20% larger than at the critical temperature 
and we may require v(T c )/T c > 1. Now, the ratio v/T at the critical temperature is a function of 
x. The use of the lattice MC results, together with perturbation theory, allows us to determine 
this function quite reliably^, and the lower limit on the vacuum expectation value is converted 
to the upper limit on the ratio of constants, x < 0.043. For an extreme opposite case, when 
A > and the Universe is reheated up to the critical temperature, the constraint is somewhat 
stronger, x < 0.026. To summarize, electroweak baryogenesis requires that the parameter X^/g^ 
in the 3d SU{2) x U(l) gauge-Higgs effective theory is bounded from above, 

X 3 /gl < 0.026- 0.043. (7.3) 

In order to obtain the constraints, following from this requirement, on the particle spectrum of 
the underlying 4d theory, one has to express this ratio through the physical parameters of the 4d 



theory at the critical temperature. This computation may be quite involved ||107|| , but it is very 
clean from the physics point of view and does not contain any infrared divergences. An essential 
point is that only one-loop graphs have to be computed in weakly coupled gauge theories, such 
as the MSM or MSSM. 

The application of the constraint of eq. (|7.3| ) to the case of the MSM follows from fig. || 
For experimental numbers m t = 175 ± 20 GeV, and win > 65 GeV one finds x > 0.07, which 
is inconsistent with eq. ( |7.3j ). Moreover, for m t = 175 GeV no Higgs mass can ensure the 
necessary requirement of eq. (I7.3D P1. So, in the MSM the baryon number non-conservation is in 
thermal equilibrium after the phase transition. This points to new physics at the electroweak 
scale, which may strengthen the first order nature of the electroweak phase transition. 

19 The complete SU{2) x U(l) gauge-Higgs model in 3d has never been simulated on the lattice, so that the 
treatment of the U(l) factor is perturbative. 

20 The upper bound on the Higgs mass was evolving in time quite a bit. The first estimate of the critical Higgs 
mass, based on the one-loop effective potential for th e sca lar field with small mass of the top quark (which was 
unknown at the time) , gave a value of mn < 45 GeV |232|, p5| . The accounting for the large top quark mass and 



Debye screening effects in the one- loop effective potential reduced this number to mn < 35 GeV [137|. The two- 



loop effects 148 |l49f somewhat relaxed this condition, while the assumption about the large non-perturbative 



effects in the symmetric phase [248] allowed a sufficiently strong first order phase transition with experimentally 



allowed Higgs boson. The lattice simulations 1 1 69 , 165| reduced all uncertainties substantially. 
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The two Higgs doublet model has more freedom, and the results of refs. IP-83] , |184fl show that 
the constraint ( [7.31 ) can be satisfied there. The reason is that the effective 3d scalar self-coupling 
constant is a complicated combination of the different scalar and pseudoscalar masses and mixing 
angles. The extensions of the standard model (supersymmetric or not), including scalar singlets, 
can also help |4|, |50|, |51 . 



According to refs. ||19l|, |193 



194 | the phase transition in the MSSM in the most part of the 
parameter space occurs in the same way as it does in the MSM. Here the MSSM also fails in 
preserving the baryon asymmetry after the phase transition. However, in a recent paper ||252| a 
specific portion of the parameter space of the MSSM, where electroweak baryogenesis is possible, 
has been founcQ What is most interesting is that quite strong constraints on the masses of the 
Higgs boson and squarks were derived. 

In order to explain the idea of ref. ||252|| in a simplest way let us add to the minimal standard 
model an SU(2) singlet but colour triplet scalar field (scalar quark) \ with the potential 



U(x, $) = —zm 2 H &$ + A($ f <l>) 2 + m 2 X *X + 2h X *X&® + ^s(x*xf 



(7.4) 



Assume now that the expectation value of the field x is zero a t a U temperatures (this is possible 
at some particular choice of parameters). Then the contribution of this field to the effective high 
temperature Higgs potential is 

-^(m 2 (T) + ^)f. (7.5) 

Now, if the effective high temperature mass m(T) is small near the electroweak phase transition, 
m 2 (T c ) ~ 0, then this term increases the magnitude of the cubic coupling a in the effective 
potential ( |5.13|) , a — > 9g 3 /(327r) + 3/12/ (27r). This, in turn, makes the phase transition stronger 
first order, and the value 4>{T C )/T C (see eq. ( p,17[ )), crucial for the electroweak baryogenesis, 
increases^. 

In the case of MSSM the role of SU(2) singlet is played by the right handed light stop |25'4 . 
Its high temperature effective mass m 2 (T) contains two essential contributions. The first one is 
the soft supersymmetry breaking mass, and the second is a positive temperature contribution 
~ g 2 T 2 , where g s is the strong gauge coupling constant. To make the idea work, the soft SUSY 
breaking mass must be negative and approximately equal to the high temperature contribution 
at the critical temperature. Previously, the negative values of that mass have not been considered 



because of the danger of colour breaking; the authors of [[252j| have shown that it is possible to 
satisfy simultaneously the requirements of the absence of colour symmetry breaking, strongly 
enough first order phase transition together with experimental bounds on SUSY particles. The 
region of parameters allowing for electroweak baryogenesis requires that the Higgs mass is smaller 
than the Z mass, lightest stop mass is smaller than the top mass, and tan (3 < 3. This range of 
masses is accessible for experimental search at LEP2 and Tevatron. 



21 We thank M. Carena, M. Quiros and C.E.M. Wagner for describing their results prior to publication. 

22 A similar idea of introducing SU(3) x SU(2) x U(l ) sca lar singlets in order to e nhance the strength of the 



electroweak phase transition has been suggested in ref. |24 
temperature contributions to the SU(2) singlet scalar mass 



A new element of [ 252 is accounting for the high 
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Yet another possibility to have strongly enough first order phase transition may be realized in 
the models of dynamical electroweak symmetry breaking [ |253| , [254j| . Here the phase transition 
occurs in a strong coupling regime both in the symmetric and Higgs phases and high temperature 
3d description does not work. These type of models predict new (strongly interacting) physics 
at the TeV scale. 



7.2 Sources of CP-violation in the EW theory and its extensions 

To produce the baryon asymmetry, the particle interactions must break C and CP symmetry. 
C symmetry is broken due to the chiral character of electroweak interactions. In the Minimal 
Standard Model, the conventional source of CP violation is that associated with Kobayashi- 
Maskawa (KM) mixing of quarks. The Yukawa interaction of quarks with the Higgs boson in 
the MSM has the following form, 

L Y = -j^—{Q L KM d D R <p + Q L M U U R 4> + h.c.}, (7.6) 

where M u and are diagonal mass matrices of up and down quarks, K is the KM mixing 
matrix, containing one CP violating phase 5cp- The MSM contains yet another source of CP- 
violation, associated with the QCD vacuum angle 9. It is constrained experimentally, 9 < 10~ 9 . 

A popular extension of the MSM is a model with two Higgs doublets, ipi and ip 2 - In order 
to suppress flavour changing neutral currents, the interaction of Higgs bosons with fermions is 
chosen in such a way that <p\ couples only to right-handed up quarks while <p 2 couples only to 
down quarks. The other possibility is that <^ 2 decouples from fermions completely and tpi gives 
masses to all the fermions. In addition to the KM mixing, this model contains CP violation in 
the Higgs sector. The scalar potential has the form ||255|| : 



V = \i((p\<pi - vlf + \ 2 [^\^ 2 - v 2 2 f+ 

A 3 [(^i - v\) + (v4^>2 - v 2 2 )} 2 + 
M{Avi){V292) ~ (A^iflfi)} + (7-7) 
\ 5 [Re(tp\(p 2 ) -ViV 2 cos£] 2 + X 6 [lm.((p\(p 2 ) - v x v 2 sin£] 2 , 

£ being a CP-violating phase. 

In the super symmetric extensions of the standard model the Higgs potential is CP invariant 
and CP is violated by the soft supersymmetry breaking terms. In the simplest version of MSSM 
there are two extra CP phases and the relevant interaction has the form | |256| , (for a review see 
ref. p57| ) 



[fiHH'] F + m g [A{UtuQH + D£ D QH' + E&LH') + {i b HH']a + h.c. (7.8) 

where U, D, Q, L, E, H and H' are the quark, lepton and Higgs superfields respectively, 
parameters /i and A are complex and flavour matrices £ are assumed to be real, m g is the gravitino 
mass. In this model extra CP-violating phases appear in the vertices containing superpartners 
of ordinary particles. 
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7.3 EW baryogenesis: how to state the problem 



Switching off the sphaleron transitions in the broken phase is, clearly, not enough for the asym- 
metry production. According to the third Sakharov condition, baryogenesis requires deviations 
from thermal equilibrium in reactions that break CP. This is provided by the first order nature 
of the phase transition, which proceeds through the bubble nucleation. Before the bubbles per- 
colate, the largest deviations from thermal equilibrium (e.g. in the particle number densities) 
are at the fronts of the shock waves of the deflagration bubbles and near bubble walls. At the 
time of percolation, deviations from thermal equilibrium arise because of collisions of the do- 
main walls and shock fronts. The latter effect cannot give substantial contribution to the baryon 
asymmetry since it is proportional to the fraction of the Universe volume occupied by domain 
walls, 

4" ~ 1(rl °> ( 7 - 9 ) 

where £ is the thickness of a domain wall (scalar correlation length), Rf, is the typical bubble size; 
for numerical estimates of £ and Rb see section 5. The shock fronts propagate in the symmetric 
phase, where the rate of fermion number non-conservation is higher than the rate of the Universe 
expansion, or in the broken phase, where B-nonconservation is switched off. So, baryogenesis 
cannot happen near shock fronts, and the only possibility is to associate it with domain walls. 

For the most part of their life, bubbles of the broken phase are macroscopic (their size is much 
larger than the typical correlation length); the domain walls move with constant velocity before 
percolation; after it the Universe may or may not be reheated up to the critical temperature, 
depending on parameters of the model. If it does, the velocity of the bubble walls drops con- 
siderably down to values v ~ 10 _2 -10~ 3 , and then slowly varies depending on the bubble size. 
Therefore, the picture of a planar domain wall, "eating up" the symmetric phase with some 
velocity v, is a good approximation to the problem. 

How should the solution of the baryogenesis problem look like? In very general terms, the 
answer is: Write down the kinetic equations accounting for all relevant processes, supply them 
with appropriate boundary conditions (equilibrium in the symmetric phase far from the domain 
wall) and then determine the baryon number deeply in the broken phase by solving these equa- 
tions. Besides the two obvious types of processes (1. B non-conservation, which is rapid in the 
symmetric phase and slow in the broken one; 2. CP-violating interactions of various particles 
with the domain wall and with each other), a number of other reactions should be taken into 
account. 

The first group of phenomena deals with B and L conserving processes. 1. Ordinary strong 
and weak interactions tend to make momentum dependence of distribution functions for quarks, 
leptons, gauge bosons and Higgses to be an equilibrium one. These processes govern the diffusion 
of the CP asymmetries in fermion number, created in the vicinity of the domain walls. 2. 
Chirality flip interactions of quarks and leptons, coming from interactions with Higgses and 
from strong sphalerons. These reactions tend to make the concentrations of left-handed and 
right-handed particles equal to each other. Since anomalous B and L non-conservation deals 
with left-handed fermions, the left-right transitions influence the 5-violating reactions. 3. Debye 
screening of the long range gauge forces, which tend to damp any non-trivial distribution of the 
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dynamical charges, such as hypercharge [|258|1 (see also ||259|| ). 

The second group of phenomena deals with the description of the relevant degrees of freedom 
at high temperatures: 1. High temperature physical excitations are different from those at zero 
temperatures. Therefore, the corresponding kinetic equation should deal with quasi-particles 
rather than particles ||147| , f260j , |124| . 2. Simultaneous interaction of quasiparticles with the heat 
bath and varying scalar field results in the mixing, analogous to that of the neutrino in matter 
26 1| , |262|| . So, quasi-particles are to be described by the density matrices rather than the particle 
number distributions [|63|, |6|, |65|, |6|. 3. Q uasi-particles in the plasma have finite lifetime, 
i.e. they should be characterized by an energy and momentum simultaneously. In general, the 
kinetic equation should be able to account for this. 

To our knowledge, the complete programme outlined above has never been carried out. The 
main difficulty is the construction of the kinetic equation incorporating all necessary features. 
Some of the effects mentioned above were taken into account, but the complete picture is still 
missing. So, we consider, at the qualitative level, various ideas and estimates of the baryon 
asymmetry produced at the EW phase transition. 



7.4 Uniform scalar fields 

A good theoretical laboratory, allowing an understanding of physical processes giving rise to 
the charge asymmetry, is the consideration of the uniform but time dependent scalar fields. 
Probably, this situation is never realized, but this case is much simpler than that of the bubble 
wall propagation. 

Suppose that we have a kind of spinodial decomposition phase transition, in which case the 
scalar field is initially near = and the system is in the symmetric phase. Sphaleron processes 
are in thermal equilibrium. Then the scalar field uniformly rolls down to the true vacuum, where 
the SU(2) x U(l) symmetry is broken and sphaleron processes are suppressed. The first rough 
estimates of the baryon asymmetry in this case were given in ref. ||267||, and a lot of work on 



this subject has been done in refs. ||268| , |269| , p70j , |271| , [272j , |273|| and many others, for relatively 



recent reviews see refs. and references therein. 

We will consider the main idea on the example of the two Higgs doublet model. Our scalar 
fields ipi and tp 2 are uniform in space but change from ip = to ip = ip c during time At of the 
spinodial decomposition phase transition. Suppose that this time is small enough, At/r top <C 1 
where T top is a typical time of top quark chirality flip (the top quark is most important since 
it has the largest Yukawa coupling constant). Then the top quark distribution has no time 
to adjust itself to the changing scalar field. So, it may be integrated out with the use of the 
equilibrium Matsubara technique. This was carried out in ref. [271] with the result that the 
effective action has the following form: 



Sp = fxN cs , (7.10) 



where 



/' - -4C(3) (^Y^Ofa), Ofa) = (<p\V <pi - (Po^i)Vi), (7-11) 



' \ixTJ vf 

and m t is the mass of the top quark, ( is the Rieman (^-function. 
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The effective bosonic action now breaks P and CP simultaneously, with CP violation in the 
scalar field potential, and P violation in the term ( |7.10| ). This allows us to generate the non- 
zero value of the topological charge Q, which is P- and CP-odd Note also that in more 
complicated models the operator /i may appear to be P- and CP-even (e.g., \i ~ d (<p^ip)). In 
the latter case the effective action ( |7.10| ) itself breaks P and CP simultaneously. Estimates of 
the baryon asymmetry produced in this case were made in refs. [ |267| . |272|| . 

The zero-temperature bosonic effective action of this model also contains parity odd term 
8(x)q(x), where 8(x) is the relative phase of the scalar fields. This important fact was discovered 
by Turok and Zadrozny and applied to baryogenesis in ref. ||2G9| . 

If At is not too small, At ■ m(T) 3> 1, where m(T) is the typical mass scale at high tempera- 
tures, then the term (|7.10|) may be considered as the chemical potential for the Chern-Simons 
number and the number density of fermions created during the transition is 



n B = n f J dtT sph (t)fi(t), 



(7.12) 



where T sp h is the time-dependent rate of the sphaleron transitions, and rif = 3 is the number 
of fermion generations. The sphaleron rate T sp h rapidly decreases when the mass of the vector 
boson increases. A natural way to estimate it in the entire range of W-boson masses is to use 
eq. (|4. 14] ) for some Mw > M crit and ([149]) for the opposite case where M crit ~ 7awT is foundQ 



from the relation T 



br 



sym- 



In this approximation, 

n B ~ 7if(a w T) 4 fx(t* 



(7.13) 



where /i(t*) is the chemical potential at the "freezing" time my/{t*) = M crit . The asymmetry was 
estimated in ref. [|271 | and lately corrected in the detailed analysis of ref. [^| for the spinodial 
decomposition phase transition. The asymmetry reads 



A 
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eff 



-KUja^r sin 3 2a\cp , 



v{v 2 



(7.14) 



with N e ff being the number of effectively massless degrees of freedom, X C p = (A 5 — A 6 ) sin2£ , 
and 

1 (7.15) 



tana 



m|(T c 



with rrii(T c ) being the temperature-dependent scalar masses at the moment of the phase tran- 
sition (see, e.g. ref. |183|). An analogous dependence on the coupling constants was found in 



ref. |[274|| . In spite of the rather high power of the coupling constants, this estimate can give an 
asymmetry consistent with observations 0. 

23 The presence of fermions is essential here. The purely bosonic tree action conserves P, and the net topological 
charge cannot be produced. 

24 A factor 2-3 instead of 7 was found in ref. |^74| from other considerations. Clearly, these estimates are 
qualitative rather than quantitative. 



25 One can obtain a similar estimate for the asymmetry from different consideration J 267 , [269 , 27C | dealing with 
non-perturbative fluctuations of the Chern-Simons number in the symmetric phase [275|. 
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This consideration can be easily generalized to more complicated models. First, one calculates 
an effective bosonic action, which breaks P and CP and defines an effective potential for the CS 
number. Then, an estimate of the net production of fermions is given by (|7.12|) . 

An essential assumption in the estimates presented above is that the time of the phase transi- 
tion is shorter than that of kinetic reactions. An "exact" solution to the problem in the opposite 
case for quite a specific situation has been suggested in ref. ||273|| . Since this example is in- 



structive, we reproduce here the main idea of this paper, using correct coefficients from ref. 
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Let us take again the two Higgs doublet model. For simplicity, we set all Yukawa couplings, 
except for that of the t-quark, to zero, i.e. the Yukawa interaction is assumed to be 

L Y = ftQsUsVi- (7.16) 

Here Qi are the left-handed fermion doublets, U^, Di are right-handed quark fields, and i is the 
generation index. We also neglect the reactions with quark chirality flip associated with strong 
sphalerons. 

We put A5 = A6 = in eq. (|7\7| ); with these couplings the potential has an extra global 
U(l) symmetry, which is spontaneously broken. Let us consider this model in an extermal Higgs 
background of a special form, namely 9 = 9t at t > and 9 = at t < 0, where 9 is the 
Goldstone mode, 

Suppose that at t < the system was in thermal equilibrium and was charge symmetric. We 
want to determine the baryon number of the system at t — > 00. The density matrix p(t) of the 
system obeys the Liouville equation 

i^ = mt),P(t)}, (7-18) 

where H(t) is the time-dependent Hamiltonian of the system in the background field. Now, 
one can make an anomaly free hypercharge rotation of the fermion fields in such a way that 
the time dependence disappears from the Yukawa coupling ( |7.16| ). Because of the global U(l) 



tanfl= ^7^ . (7.17) 



symmetry this converts the time dependent Hamiltonian to a time independent one, H{t) 
H e ff = H — 9Yp where Yp is the fermionic hypercharge operator, 



1- 4 - 2 - 

-^QiloQi + -UijoUi - -A70A - Li<y Li - lE^Ei 



(7.19) 



At t — > 00 the system must be in thermal equilibrium, = 0. Since in the new representation 
the Hamiltonian is time independent, the density matrix is 



p(oo) = - exp 



(7.20) 



where Xi is a complete set of conserved charges (operators commuting with the Hamiltonian). 
Their average must be equal to zero. This requirement fixes the chemical potentials /i, and 
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allows the unambiguous determination of the baryon number of the system. A complete list 
of the conserved charges can be found in ref. [ |125|| , and we quote here the final result for the 
baryon number only, 

(B) = ^_I*6>(1 + OK 2 /t 2 )), (7.21) 
where n s = 2 is the number of the scalar doublets. 



The result ( [7.2 1|) is quite amazing. It does not contain the Yukawa coupling constant, the 
scalar vacuum expectation value, or the rate of sphaleron transitions. One might even think that 
it is wrong, since if f t or v or T sp h is zero, then, obviously, one must have (B) = 0. Nevertheless, 
it is correct. The key point is that the time at which the asymptotic value of the baryon number 
is reached tends to infinity when the above mentioned quantities tend to zero. Many conclusions 
based on the straightforward analysis of the perturbation theory break down at large times, 
when the application of the kinetic theory is essential, and this is one of the examples. For 
typical values of the parameters, the top quark chirality equilibration time is r t ~ 30/T, and 
the B non- conservation time is T sp h ~ 10 5 /T; the result ( |7.21| ) is valid only for t > T sp h- The 



discussion of intermediate time r t <C t <C r sph is contained in ref. |[273| , [L25| 



It is worth noting that high temperature sphalerons, and other chirality flip reactions, may 
change the estimate ( [7.21| ). In particular, strong sphalerons, discussed in Section 4 have the 
physical effect of maintaining the same chemical potential for left- and right-handed baryonic 
numbers and diminishing the set of conserved quantum numbers in the system. This leads to 
the suppression of the baryon number by a factor ~ (^) 2 ||125|| . Other aspects of influence of 



strong sphalerons on baryon asymmetry were discussed in ref. 

We conclude this discussion by remarking that the use of U(l) global symmetry was essential 
in the derivation of eq. ( |7.21| ). Without it the hypercharge rotation would not remove the time 
dependence from the Hamiltonian, and the solution of the Liouville equation could not be found 
so straightforwardly. This is discussed in more detail in refs. [|274|, E77I . 



7.5 Asymmetry from fermion— domain wall interactions 

In reality, though, the phase transition goes through the bubble nucleation rather than as spin- 
odial decomposition. This is an additional challenge, since the baryon number (or, in general, 
asymmetries in particle number densities) can now be distributed in a non-uniform way and 
depend on the distance from the domain wall. Correspondingly, the analysis of the kinetic 
equations is much more complicated. 

Two different cases are usually considered, depending on the relation between the mean free 
paths of particles and domain wall thickness. The physics of thick wall baryogenesis was orig- 
inally considered in refs. [271 , [272j| and has much in common with the quasi-adiabatic case 



of uniform fields discussed in the previous subsection. P or CP non-invariant interaction of 
fermions with the moving domain wall together with CP-breaking scalar dynamics induces P- 
and CP-odd terms in the bosonic effective action, which bias the sphaleron transitions inside 
the domain wall. The excess of quarks generated in this way is absorbed then by the expanding 
bubble. (Another, equivalent way to say this [ 275 ] is that the particle densities of fermions 



gradually adapt to scalar background changing in space and time in such a way that an excess 
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of right quarks and left antiquarks is created. Left antiquarks are destroyed by the sphaleron 
reactions while the right fermion number is intact and is converted into baryon number at the 
end.) 

A nice physical picture of thin domain wall baryogenesis was suggested by Cohen, Kaplan and 
Nelson [ |278| , [279| , |280|| . Since the masses of fermions are different in the symmetric and broken 



phases, they scatter on domain walls (are reflected or transmitted). CP violation manifests itself 
in different reflection coefficients for particles and antiparticles. So, the moving domain wall 
acts like a separator for different types of fermion numbers, filling the bubble with fermions and 
outer space with antifermions (or vice versa, depending on the sign of CP violation). Of course, 
the interactions of fermions with domain wall conserve the fermion number, i.e. the number of 
fermions flying into the broken phase is equal to the number of anti-fermions moving into the 
symmetric phase. Antifermions, injected into the symmetric phase, participate in the anomalous 
reactions that change fermion number, while fermions injected into the broken phase do not. As 
a result, non-zero baryon and lepton asymmetries are established in the broken phase. 

Clearly, reliable calculations of the effect in realistic theories are quite complicated because of 
the large number of different particle species participating in interactions. Moreover, a number 
of effects, discussed above, should be taken into account. A number of papers is devoted to 
the study of the origin of the CP-violating fermion currents and their propagation in front of 



domain walls ||28 1| , [282| , |283| , |284|| ; the most recent (and probably most elaborate) treatment can 
be found in [f285|, |286|, pTTj pg7|, M- 

Below we discuss the qualitative features of the domain wall baryogenesis. Our consideration 
is by no means complete, and the reader may consult the original papers for details. 

At sufficiently small velocities of the domain walls we can divide the problem in two parts 0. 
The first one is the microscopic calculation of various fermion currents at the domain wall. The 
second one is the consideration of the diffusion of the particle number densities in front of the 
wall and their dissipation in different processes. 



We begin with the first part [[278] , |279| . [280|| . Let us ignore for a moment any high temperature 



effects. The simplest case is that of the thin domain wall moving with some constant velocity v. 
Let us choose the rest frame of the wall and consider scattering of fermions on it. For example, 
left-handed fermions incident in the symmetric phase may be reflected back to the symmetric 
phase as right-handed fermions (because of the spin conservation) or can be transmitted through. 
The transmission and reflection coefficients (i is the label of an incident fermion and j is that 
of the final state) can be found from the Dirac equation: 

with appropriate boundary conditions^]. Here L and R correspond to up and down components 
of two-dimensional Weyl spinors. The x dependent matrix M is complex, giving rise to CP- 
violation. In general, r^- for particles are different from f%j for anti-particles. This leads to 

26 If the diffusion tails (see below) are comparable with the thickness of the domain walls, this is not possible. 
27 A method for high precision numerical evaluation of reflection coefficients was considered in ref. [f24 . 
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non-zero fermionic currents: 



= f ^A(u f (uj + ) - n F (cu„)) [rV - ftf] . , (7.23) 

where n l F is the Fermi distribution for the incident particles, u± = uj ± vp t , p t and p\\ are the 
momenta of fermions tangential and parallel, respectively, to the domain wall. This expression 
vanishes if the domain wall is at rest (v = 0) or if there is no CP violation. The total baryonic 
current Jqp = Jl + Jr which results from the solution of the Dirac equation (|7.22|) vanishes, 



but the currents of left-handed ( Jl) and right-handed ( Jr) fermions are non-zero. 

The construction of the Dirac equation for quasiparticles, accounting for leading high temper- 
ature effects, was done in ref. ||260| , |124|| , where generalized expressions for the particle currents 



can be found. The major qualitative effect of high temperature corrections is that the currents 
of left-handed and right-handed fermions do not compensate each other and the total baryon 
current is produced. Physically, this happens because left-handed particles participate in the 
weak interactions but right-handed particles do not, and Jqp ~ cx-w^l- 

The thin wall description of the fermion scattering is applicable only if the mean free path of 
fermions at high temperatures is much larger than the domain wall thickness. This allows us 
to use distribution functions of fermions undisturbed by the domain wall and impose ordinary 
boundary conditions for the scattering problem at spatial infinity. The thick wall case (the 
mean free path is small compared with the thickness) is much more complicated. Clearly, the 
scattering description is not adequate in that case. The physical phenomenon to be taken into 
account is the modification of the particle distributions across the domain wall. A number of 
interesting effects arising in the latter situation are discussed in ref. [[286| , |277| , |285|| . 

The second problem is the particle transport. Several approaches were applied to the consid- 
eration of it. The first one is that of Monte Carlo simulations of the injected flux of particles 
280|| , the second one is diffusion equations |12|, |77|, |287|, |8|, [189|, |288|. Limitations of the 



diffusion approximation were considered in ref. [p89j. In the discussion below we closely follow 



ref. ||124|| where the analytical approximation to the problem was constructed for a simple case. 

Given the flavour and chirality structure of the fermionic currents, the diffusion equations 
should be written for all particle species. Fermions participate in many processes on both sides 
of the wall with different time scales. In order to understand what the relevant time scales are, 
let us consider the fate of a particle after it has been reflected from the domain wall towards 
the unbroken phase. Roughly, its typical distance from the bubble wall is given by \f~Dt — vt, D 
is the diffusion coefficient. The first term describes the random walk of the particle in the rest 
frame of the plasma and the second term describes the motion of the bubble wall. This particle 
will be trapped by the bubble after the time interval to ~ D/v 2 , so that all processes with 
characteristic time r < must be taken into account. The examples of the relevant processes 
include 5-violation, the elastic scattering of quarks and gluons, chirality flipping transitions of 
heavy quarks, strong sphalerons. 

In order to get a better feeling of the physics involved, consider the simplest case when the 
total baryonic current originated from CP non-invariant interactions is not zero and neglect all 
processes besides the elastic scattering of fermions and anomalous B and L non-conservation. 
Let us take a planar domain wall which moves through the plasma with sufficiently small velocity 
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v (we shall see below how small it should be in order that this consideration works). We take 
a reference frame associated with the domain wall; let the broken phase be to the right and the 
symmetric one be to the left, and x be the distance from the domain wall. Assume that the 
thickness of the domain wall is small enough (again, we shall see below what this means). We 
denote by n B (x, t) and n^x^t) the densities of baryon and lepton numbers in the rest frame of 
the wall. The diffusion equations for the broken phase, where sphalerons do not operate, are 






(7.24) 



For x < we have 




D B £^-v4z-lv -r 



9 2 _ ,,d_ _ 3 
U dx 2 " 




(7.25) 

where T = 9T sph /T 3 , D B and D L are diffusion constants for quarks and leptons respectively. An 



estimate of these gives ||277| , |287fl D B ~ Di ~ 



We are looking for a steady state (time independent) solution to these equations. In the 
broken phase the only solution consistent with the boundary conditions is constant density, 

n B = n L = const = B + . (7.26) 

In the symmetric phase, the solution is a combination of dying exponentials. We present it in 
two limiting cases. The first one deals with "large" velocities, 

p = 3D B T/v 2 < 1. 

Physically, this corresponds to the situation when an extra antibaryon, injected into the sym- 
metric phase, is exposed to the S-violating reactions for short times t ~ D B /v 2 , before it is 



trapped by the moving domain wall. One finds ||124|| 



(vx \ f vx \ 

~Db~) ' nL = C ' 2exp \~dI) ' ( 7 ~ 27 ^ 

corresponding to the diffusion of quarks and leptons to the distances D B /v and D B /v respec- 
tively. In the opposite case of low velocity, the transitions from quarks to leptons due to sphaleron 
processes are essential, and the solution reads 



n B = C 3 exp (||-) + exp ^JjL*) , (7.28) 
3 (vx\ D B 




One of the requirements for the validity of the diffusion approximation is that the diffusion tail 
(the shortest one is for the quarks) is much longer than the domain wall thickness I, namely 
I < D B /v. 
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The constants C1—C4 can be determined from the matching conditions at the domain wall. 
If we denote by Jcp the total CP-odd baryonic current originating from interactions of quarks 
with the domain wall and assume that CP asymmetry in the leptonic current is zero, then 

12 

B + = —J C pfs P h{p)i (7.29) 

where f sph (p) = 1 for p > 1 and f sph {p) = \p for p < 1. 

The asymmetry inside the bubble has an interesting velocity dependence. If Jqp ~ v, as in 
the quantum-mechanical consideration of the thin wall case, then the maximum asymmetry is 
produced at p ~ 1, i.e. v ~ y/^YDs — 0.01. It is worth noting that these small velocities 
are quite possible at the final stage of the phase transition if the Universe is reheated up to 
the critical temperature. The analysis of the consequences of this scenario can be found in ref. 
[ 229|| . For a more realistic thick wall case one effectively has Jcp ~ v 2 ||286| , [277] , |287|| 0, and the 



asymmetry is velocity independent. The same conclusion has been reached also in ref. 

In the realistic case of many particle species this consideration must be generalized. Instead of 
the CP-violating flavour independent baryonic current considered above, many CP-odd currents 
appear, resembling the flavour dependent interaction of fermions with the domain wall. The 
left- and right-handed currents must be distinguished, since particles of different chiralities have 
different interactions with the heat bath and sphalerons. Quantitatively, the results are model 
dependent. For example, in some schemes the lepton interactions with domain walls produce 
more asymmetry than quark interactions ||290|| . Serious investigations of realistic models have 



been carried out in very interesting papers [ [28 5| , [287| , [277| , |189|| , and we refer to them for more 
detail. 



7.6 Strength of CP violation and baryon asymmetry 

Extensions of the standard model, having strong enough first order phase transition, may not 
contain any new source of CP violation. The question arises whether the KM source of CP 
violation (or QCD vacuum angle) alone can be responsible for the baryon asymmetry. Let us 
begin with the KM mechanism of CP violation^} 

An important property of the interaction ( |7.6|) is that the CP- violating phase 5c p can be 
rotated away by phase transformation of the fermion fields if there is a degeneracy in the up 
or down quark sectors, or if some mixing angle between different generations is zero. In other 
words, in the standard model, CP violation vanishes together with the Jarlskog determinant 
29Tj , 

d C p = sin(6 ) i 2 ) sin(# 2 3) sm(6>i 3 ) sin 5 C p* (7.30) 

/ 2 2\t 2 2\f 2 2\ I 2 2\i 2 2\ i 2 2\ 

[m t - m c ){m t - m u )(m c - mj ■ [m h - m s )(m b - m d ){m s - m d ), 
where 8^ are the mixing angles and are the quark masses. 
28 We thank Michael Joyce for clarification of this point. 

29 Of course, this possibility is rather unnatural: the only known way to strengthen the phase transition is to 
add extra scalar particles. This means that new scalar interactions appear. In general, they contain complex 
phases and lead to CP violation. 
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The structure of the KM mechanism of CP violation makes baryogenesis a very non-trivial 
problem. Indeed, the electroweak phase transition, where strong deviations from thermal equi- 
librium are expected, occurs at temperatures of order 100 GeV. It seems, therefore, that quark 
masses (maybe with the exception of the t-quark) can be treated as perturbations, so that the 
dimensionless measure of CP violation is just 5qp ~ d C p/T 12 ~ 10~ 20 [[23 1| , [232|| . Clearly, this 
number is too small to account for the observed asymmetry However, there may be loopholes in 
this argument and a number of dynamical mechanisms, in which the KM source of CP- violation 



may be enhanced, have been suggested [|232| , |260| , |1 24 



If there is a dynamical spontaneous CP violation in the electroweak theory before [232 



or 



during [|292| , p84j] the electroweak phase transition, then the Universe contains domains with 



different CP-parity at some stage in its evolution. The small explicit CP violation breaks the 
degeneracy between the different CP states, so that more energetically favorable domains "eat" 
those with the opposite CP parity. The fact that the age of the Universe at T ~ 100 GeV is 
macroscopic, gives rise to an enhancement factor ~ Mpi/T* ~ 10 16 [232 . In this mechanism, 



the baryon asymmetry does not depend on the magnitude of CP violation but does depend on 
and a power-counting estimate of the effect gives 



its sign 



292 



A ~ — a 

Neff 



3 

W ■ 



(7.31) 



We should stress, however, that the possibility of spontaneous CP breaking at high temperatures 
in electroweak theory is very speculative. 

Another possible caveat in the no-go argument presented above has been discussed in ref. 
124j| . The "no-go" theorem relies on the applicability of the perturbation theory in quark 
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masses and makes use of the assumption that the typical energy scale relevant to the estimate 
of the asymmetry is the temperature of the phase transition. In fact, these assumptions break 
down when the interaction of fermions with domain walls is considered. Namely, if the energy 
of the quark in the unbroken phase is smaller than its mass in the broken phase, it will be 
reflected from the domain wall with unit probability independently of the value of its mass. The 
perturbation theory does not work only in a small fraction of the phase space determined by the 
quark masses, but the loss in the phase space factor may be smaller than the gain in CP-violating 
amplitude. The maximum asymmetry comes from the strange quark complete reflection: the 
GIM cancellation does not occur for it, and what is left over from dcp is just the product of 
mixing angles and CP- violating phaseQof order 10~ 5 . The estimates of the asymmetry presented 
in ref. ||124j| are rather uncertain, A ~ 10 -10 -10~ 18 , but they indicate that the KM mechanism 



of CP violation cannot be discounted as a source of the baryon asymmetry of the Universe. 

The possible impact of the strong CP violation (associated with QCD vacuum angle 9) on the 
electroweak baryogenesis has been analysed in ref. ||299|| . It was concluded that it does not play 
any significant role because strong CP effects at the electroweak phase transition temperature are 
suppressed by at least semi-classical exponent exp(— 2p) ~ 10~ 25 , yet other suppression factors 
come from Yukawa couplings. The same statement holds true in the models with axions, where 
effective angle 9 may be of the order of 1 at the electroweak scale. 



30 In refs. [293, 294, 295, 296 1 an opposite conclusion was reached. The calculational procedure of these works 
was criticized in ref. [297, 298|, where it was argued that the claims of ref. [293, 294, 295, 296 1 are not justified. 
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While convincing arguments showing that the standard model CP violation is sufficient to 
generate the baryon asymmetry are absent, the extensions of the electroweak theory naturally 
provide new sources of CP violation. In spite of the fact that the existing estimates of the baryon 
asymmetry are valid, probably, within an order of magnitude, it is clear that extended versions 
of the electroweak theory can accommodate the observed baryon asymmetry of the Universe. 
The specific estimates for the two Higgs doublet model can be found in ref. [ 285 , [277] , 287 . 189] 
and for the sypersymmetric theories in ref. 



300, 301, 302, 303 



8 Instanton-like processes in high energy collisions. 

As discussed in section 2, instanton-like transitions may occur at unsuppressed rates at suffi- 
ciently high energies. This possibility is definitely realized at high temperatures. It is natural 
to ask whether collisions of highly energetic particles may lead to baryon and lepton number 
non-conservation with exponentially unsuppressed cross sections. The relevant energy scale, 
E sp h ~ ray/ /aw (10 TeV in the electroweak theory), is not too far from collider energies, so this 
problem is not of academic interest only. 

In this section we outline the current status of this very complicated problem. We summarize 
the results of perturbative calculations about the instanton, which show that the cross section 
indeed increases exponentially with energy at E <C E sp h. However, the perturbation theory 
about the instanton is unreliable in the most interesting energy region E ~ Esph, so the pertur- 
bative calculations cannot tell whether the instanton factor ( |2.8| ) is overcome. We present below 
a rather general argument based on unitarity and conventional perturbation theory, which shows 
that the instanton-like transition rates are most likely exponentially small at all energies. This 
argument, however, does not exclude the possibility that the suppression disappears asymptot- 
ically at E — > oo, and that the actual values of the suppression factor are not too small at 
realistic values of coupling constants and energies. So, the computation of the instanton-like 
transition probabilities remains an interesting problem whose solution requires non-perturbative 
approaches. We outline in this section one of these approaches and first non-perturbative results; 
these results indicate that the instanton-like cross sections are indeed unobservably small at all 
energies. 

In spite of the considerable progress in understanding the instanton-like processes in high 
energy collisions, this problem is still not solved completely. The non-perturbative techniques 
adequate to this problem are still being developed. It is worth pointing out that electroweak 
B and L non-conservation in high energy collisions belongs to a wider class of processes, which 
includes false vacuum decay induced by particle collisions, induced decays of metastable solitons 
||304j| , and, notably, production of multiparticle final states in the trivial vacuum sector (without 
instantons) [j305| , [306, 307| , 308 , 309|| . We shall not discuss the latter problem, which is also of 



potential phenomenological interest, and refer to a review ||310|| . However, we stress that its 
solution also requires novel non-perturbative techniques which should have much in common 
with the approaches relevant to instanton-like processes. 
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8.1 Summary of perturbative analysis about the instanton 

The perturbative analysis about the instanton was started in refs. |39| , [40(1 and has lead to the 
picture of exponentially increasing total instanton-like cross sections at relatively low center- 
of-mass energies, E -C E sph . It also suggested the functional form of the cross section which 



indicated that the cross section may be calculable in a semiclassical way ||311| , |312| , |313|| . We 
do not consider the technical details of the perturbative calculations which are reviewed in refs. 
314 , |315| , pi6|| and present only basic ideas and main results. 



Consider a process where two M^-bosons scatter into n V^-bosons, and the system simultane- 
ously makes the transition from one vacuum of fig. 1 to a neighbouring vacuum. The topological 
number N[A] of the relevant field configurations, given by eq. (|2.3|) , should be equal to 1 (in 
fact, in the case of a finite number of incoming and outgoing particles it is more appropriate to 



measure the change in the vacuum number n[u] by the winding number of the Higgs field ||317|| ). 
Let us disregard fermions — they play a minor role as far as the cross sections are concerned 
||318|| — and consider the bosonic sector of SU (2) theory with one Higgs doublet, i.e. the bosonic 



sector of the simplified model introduced in section 2. To evaluate the amplitude of this process 
we begin with the (2 + n)-point Euclidean Green's function 

G n+2 (x 1 ,x 2 ,y 1 ,...,y n ) = / DA e- s[A] A(x 1 )A(x 2 )A(y 1 ) . . . A{y n ) (8.1) 



where spatial and group indices are omitted. Since the instanton is a minimum of the Euclidean 
action (in fact, we have to deal with constrained instantons considered in section 2), we make 
use of the semiclassical approximation and obtain in the leading order 

Gftfix!, . . . , y n ) = [ d'x % ^{p)e-^ 2p2v2 A ms \ Xl -x ;p)... A ms \y n - x ; p). (8.2) 

J p 

Here xq and p are the usual collective coordinates of the instanton and the measure is the same 
as in eq. fl2.9p. The constrained instanton configuration is described by eqs. ( J2.1UQ and fl2.7|). 
One should also integrate over the instanton orientations; this integral is not explicitly written 
in eq. (|8.2|) . The Green function ( |3.2|) may be picturized as shown in fig|| It is worth pointing 
out that the dependence on coordinates in eq. ( |8.2| ) factorizes, up to the integration over the 
instanton position ensuring the overall momentum conservation. This means that, in the leading 
semiclassical order, the Green function has a point-like structure. So, in this order, the cross 
section of the process 2W — > nW will exhibit a power law increase with energy. 

To obtain the cross section, one performs the analytical continuation into Minkowski space- 
time and then makes use of the LSZ procedure. This can be done easily, as the only dependence 
on coordinates in eq.( |S.2D is through A mst . Since the constrained instanton field decays in 
Euclidean space-time according to eq.( |2.7| ), its Fourier transform indeed has a pole at p 2 = — m^. 
The residue at this pole can in fact be obtained directly from eq. ( p,10| ): the Higgs mechanism 

l w 

up to small corrections. In this way one obtains the residue 



shifts the pole of the instanton field from p 2 = to p 2 = —vnL, but does not change the residue, 



R(p;p) = VlPl (8-3) 
9 
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where we omitted the tensor structure depending on the instanton orientation. Thus, the am- 
plitude has the following form, 



i4K v (k 1 ,k B ,p 1 ,...,p fl ) 



dp _ 8 "" 2 _ 7r 2„2 1) 2 

— fjL{p)e ~^ R(k 1 -p)...R(p n ;p)8(k 1 + k 2 -pi-...-pn) (8.4) 

P 5 

where ki, k 2 and pi, . . . ,p n are the momenta of incoming and outgoing particles, respectively. 
The integration over p is straightforward, while the integration over orientations, implicit in eq. 
(fS.4|), is quite complicated. Ignoring the latter complication we obtain the following estimate for 
the amplitude, 



8 ^2 / 1 

A l 2 Zt n9 ~ e~~r(n + 2)! - |ki||k 2 || Pl | . . . \p n \5(K- Pl - . . . - p n ) (8.5) 



where K = (E,0) is the total centre-of-mass four- momentum. Equation (|8.5|) leads to the 
following estimate for the 2W —>■ nW instanton-like cross section at n ^> 1, 



a. 



I tin ling (J^\ ^ f TT ^ ^* |^2eod*Ti5|2 

(fci • fc 2 ) ™ ! H 2^ Pi (27t) 3 2 ^ n 



1 / const • E 2 



2n 



e ~a r 



n\ \ gv 2 n 

where we assumed for simplicity that the outgoing particles are relativistic, i.e. E/n ^> m^. As 



expected, the cross section exhibits the power law growth with energy |39], . The total cross 
section, in the leading semiclassical order, grows exponentially, 

leading i 7-1 \ leading 



a l 2 ZZ^E) = T,^ZT 9 (E) 



oc exp 



16vr 2 / E 4 



+ const 



g 2 \5 ,2 i' 4 



(8.6) 



where the number of particles at which this sum is saturated is of order 



1 



The exponential growth of the cross section was found first ||319|| in the context of multi-Higgs 
final states, which have turned out to be subdominant at relatively low energies, where the leading 
order calculations are reliable. This type of behaviour is inherent in all models with instantons 
||315|| and also for the processes of multiparticle production without instantons ||320| , |321| , |310| . 



In the case of 2W — > nW processes, there exist at least two ways of actually calculating the 
cross section including the constant in eq. (|S.6|) . One is to make use of instanton-anti-instanton 



61 



configurations [ 322 , [323|1 and another is based on the coherent state formalism [|3 1 If . Both 
techniques lead to the same result: 




leadinq i 
02^any{E) OC eX P 

where Eq = \fQnmw / aw ~ 15 TeV is of order of the sphaleron energy in the electroweak theory. 
In writing eq. (|8.8|) we made use of the fact that 



1 ( 1 



4 I 9 A f 

a w E$ 

For the same reason, the number of W^-bosons produced can be written as follows, 

1 / E \t 

n — . (8.9) 

aw \Eq/ 

Note that the average energy per outgoing particle is of order 

E [E\-^ 

p ~m w [ — . (8.10) 

n \E J 

Therefore, at E -C Eq the final particles are relativistic, while at E ~ Eq they are soft, E/n ~ 
raw- Note also that the typical instanton size is of order 

Jn 1 / E \ i , 
P~*Er~ — (-Er) ■ ( 8 - n ) 



E m w \E 



This estimate follows from eqs. flSTq ), ( gg ) and (JSty. 

We conclude that the leading order total cross section becomes unsuppressed exponentially at 
E ~ E'o, and at these energies the number of final particles becomes of order 1/aw- 

Clearly, the actual cross section of the instanton induced process 2W — > nW should not be 
described at all energies by the leading order result (|8.8|) : at E > Eq this expression contradicts 
unitarity. So, corrections to the leading order formula must be large at least at E > E . These 
corrections appear when the gauge field in eq. (|8.1|) is written as 



A = A inst + 5 A 



is 



and 5A is treated as quantum field. In the next-to-leading order, the action in eq. 
quadratic in 8 A, n fields in the integrand remain A mst and two fields are 5 A. Upon integration 
over 5 A one obtains the first correction to the Green function, which is similar to eq. (|8.2j) 
but with two of the instanton fields substituted by the propagator in the instanton background, 
Di ns t(z — xq,z' — x ) where z, z' are any two of the coordinates x\, X2, yi, ■ ■ ■ , y n - The first 
correction to the amplitude is then determined by the residue of the Fourier transform Di nst (q, q') 
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-m 2 W) q' 2 



-rriyy The diagrammatic representation of this correction 



at the double pole q 2 
is shown in fig. [|. 

There are basically three types of corrections, which we discuss in turn. 

i) Soft-soft corrections. These appear, in the next-to-leading order, when two instanton fields 
in eq. (|8.2|) corresponding to final particles are substituted by the propagator, as shown in fig. 
|S]a. As compared with the leading term ( |8.4| ) this contribution is suppressed by g 2 because A mst 
is proportional to 1/g while Di nst is 0(g°), but is enhanced by the combinatorial factor n 2 /2, 
the number of legs in fig.^]a that can be joined. The residue of the propagator at the double pole 
at low momenta is of order ||324| , |325| , |326jl 

,2 
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p 2 p 2 



R(p i ;p)R(p j ;p) 



where R is the residue of the instanton field, eq. (|8.3| ). Combining all factors and recalling eqs. 
( |S.9|) , ( |3.10|) and ( |3.11|) one finds that the soft-soft correction to the amplitude at relevant n is 
of order 

1 / E 

pleading I 

^Eq, 
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(8.12) 



We see that this correction exceeds the leading order amplitude even at E <C Eq. However, it has 



been shown ||311| , [312| , |313|| , that the soft-soft corrections to the total cross section exponentiate, 
so that the total cross section with these corrections included has the form 



leading 
*^2— > any 



soft— soft 



oc exp 
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(8.13) 



where we inserted numerical a coefficient 9/16 calculated in refs.| |327| , |324| , |325| , |326 | 

Higher order soft-soft corrections have also been shown to exponentiate so the total cross 
section with all soft-soft corrections included has the following functional form [ 311 , 312 , 



leadinq + soft— soft 
^any OC exp 
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(8.14) 



where E(E/E Q ) is an unknown function which, at small E/E , is represented by a series whose 
first terms are given by eq. ( |S.13|) . It has been found ||311|| that the soft-soft contributions into 
the exponent F(E/Eq) come from tree diagrams about the instanton, while loops contribute to 
the pre-exponential factor only. 

ii) Hard-hard corrections [|328| , 329 1. They are due to diagrams of fig. ^p. There is no 
combinatorial enhancement of these diagrams, but they produce large contributions nevertheless. 
The reason is that the residue of the propagator is large at high momenta |k| ~ E, 

g 2 p 2 (h ■ k 2 ) \n(h ■ k 2 )R(k i; p)R(k 2 ; p). 
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Therefore, the first hard-hard correction into the amplitude at relevant n is of order 

1 /£\Lo, „o w„„^„ 1 



j^hard— hard ^ 
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This property holds for a wide class of models 1 33C ] 
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which is again large even at E <C Eq. Higher order hard-hard corrections come from loop 
diagrams like those shown in fig. [II]. There exist strong arguments showing that hard-hard and 
hard-soft corrections exponentiate ||329| , [331|| , i.e., the total cross section has the form (|8.14|) . 



Unfortunately, the complete proof of the exponentiation is still lacking. 

iii) Hard-soft corrections. These come from diagrams of fig. |9]c. They contain both the 
combinatorial and energy factors. The estimate analogous to eqs. ( |8.12| ) and ( |8.15| ) is, up to 
logarithms, 

8 

1 / E \ 3 
j^hard—soft ^ pleading I \ 

aw \E J 

As we already pointed out, these corrections are also likely to exponentiate. 

To summarize, the perturbation theory about the instanton strongly suggests that the total 
cross section has the exponential form 



4ZL y oc exp 
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(8.16) 



The exponent is perturbatively calculable at E/E <C 1, where it is represented by a series in 
(E/E ) 2/3 (up to logarithms), 

„ / E \ 9 / E \ I 9 / E \ i . „, 

Hard-soft corrections contribute at the order (E /E ) 8 ^ 3 ; this order has been studied in refs. 
33 1 , |332| , |333|| . Hard-hard corrections begin at the order {E / E ) 10 ^ 3 . While soft -soft corrections 



to the exponent F come from tree diagrams about the instanton, hard-hard and hard-soft 
contributions include all loops. 

Clearly, the series flS.17|) blows up at E/E ~ 1. Therefore, the perturbative calculations 
about the instanton cannot tell whether the exponential suppression disappears at E ~ Eq 
or if it persists at all energies. The analysis of this most interesting problem requires entirely 
non-perturbative techniques. 



8.2 Unitarity 



Before presenting non-perturbative approaches to the evaluation of the cross sections of instanton- 
like processes, let us give a general argument in favour of the exponential suppression of these 
cross sections at all energies except, maybe, exponentially high ones. This argument is in the 



spirit of ref . ||334|| and is based on unitarity and conventional perturbation theory at low momenta 
(see also refs. [ |335| , |336| , |316|| ). 



Let us consider the full propagator of the jy-boson as shown in fig. ITT], and take the momentum 



of virtual W to be Euclidean and small, say Q 



m 



Then the dispersion relation relates this 



propagator to the total cross section of "z/ e e-annihilation" at cm. energy yfs into an arbitrary 
number of H^-bosons, 

O-tot(s) 



G{Q 



~ m 



w 



ds 



s + Q 5 



(8.18) 
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Here "z/ e e-annihilation" means just the production of a virtual VT-boson by an external probe, 
as shown in fig. [12. 

The left-hand side of eq. ([3.18|) is believed to be a nice asymptotic series in aw whose 



finite number of terms, say, k <C 1/aw, are given by ordinary perturbation theory (perturbative 
diagrams in fig. (TT|) . Indeed, at low Q 2 ~ there is no reason to suspect that non-perturbative 
contributions like instanton-anti-instanton shown in fig. |TT| are not exponentially suppressed. 
Thus, 

G(Q 2 ~ m 2 w ) =t^4 + 

i=l 

where C\~ are determined by conventional perturbation theory and k is fixed in the limit aw — ► 0. 
These k terms are precisely matched by perturbative contributions to the total cross section on 
the right-hand side of eq. (|8.18|) (perturbative graphs in fig. which include the production 
of k final particles or less. So, the instanton contribution into the right-hand side of eq. ( |8.18| ) 
is small, 

(j inst ( s ) 1 

ds tot \' < const • a k w at Q 2 ~ wL, k < . (8.19) 

s + Q l ' a w 

This certainly excludes the possibility that the instanton-like cross sections are large (of order 
a w with finite n) at energies of order E ~ mw/aw- It is straightforward to generalize this 
argument to collisions of two real vector bosons and other processes with a small number of 
incoming particles. In all cases the relation like eq.( |8.19| ) must hold. 

This argument is consistent with the expected functional form of the total instanton-like cross 
section, eq. (|8.16|) . It shows that the exponent F is negative at all energies, i.e. the instanton- 
like processes are always exponentially suppressed. This general argument does not, however, 
exclude a still very interesting possibility that F tends to zero as E — > oo, in which case the 
cross section may not be numerically small at still reasonable energies. In any case, a theoretical 
understanding of an actual exponential behaviour is definitely of interest. 

Let us note in passing that the same argument implies the exponential suppression of the 
production, in a trivial vacuum, of a large number of final particles, n ~ 1/g 2 . These processes 
have been actively studied in recent years; for a review see ref. [ |310| . 

Except unitarity, the above argument assumes the validity of the ordinary perturbation theory 
as an asymptotic expansion in aw for few-point exact Green's functions at low Q 2 . So, there 
remains a logical possibility that the perturbation theory is badly wrong starting at some finite 
order in aw- We shall see below that this logical possibility, which would be revolutionary for 
the entire quantum field theory, is not supported by existing (albeit limited) calculations. 



8.3 From many — ► many to few —> many 

The exponential form of the instanton-like cross section looks semiclassical, the failure of the 
naive semiclassical procedure being reflected by the fact that the inverse coupling constant enters 
not only the overall factor 4n/aw but also the characteristic energy scale E ~ mw/otw- So, 
it is natural to expect that there exists an appropriate modification of the naive semiclassical 
procedure, which would enable one to evaluate the exponent of the cross section. Since F(E/E ) 
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is determined both by tree graphs about the instanton (soft-soft contributions) and by loops 
(hard-hard and hard-soft contributions), the correct semiclassical procedure should incorporate 
the relevant part of loops. Clearly, the very existence of a semiclassical scheme, which would 
incorporate loops is far from obvious, and it is a challenging problem to invent such a scheme. 

The latter difficulty may be rephrased in the following way. At energies of interest, the final 
particles are soft — their average energy per particle does not contain the large parameter 1/aw 
- and numerous, n ~ 1/aw- So, it is natural that these particles may be described in classical 
terms: roughly speaking, these are classical waves. On the other hand, the initial particles are 
few in number, and carry a momentum proportional to the large parameter 1/aw- The problem 
is that these energetic initial particles are hard to describe in (semi) classical terms. 

On the basis of the above observations, the following approach to few — > many instanton- 
like transitions was suggested 



33 7 , p38|| . As an intermediate step, consider many 



many 



transitions where the number of incoming particles is 
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where v is a variable parameter. In the limit 



aw 



0. 



= fixed, 

£/ 



v 



fixed, 



5.20) 



the number of incoming particles is large, their energy per particle is of order mw /v — inde- 
pendent of aw, so these particles, as well as outgoing ones, may be described in (semi)classical 
terms. We shall see that the total probability for the optimium choice of the initial state at 
given Hi has the exponential form 



a (E,rii) oc exp 
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(8.21) 



and the exponent F(E/Eq, v) is semiclassically calculable. 

One argues that F{E / Eq,v) decreases (becomes more negative) when v (i.e. the number of 
incoming particles) decreases. Indeed, processes with fewer incoming particles may be viewed 
as a subset of processes with larger n% with some incoming particles not participating in the 
scattering. So, many — > many transitions provide an upper bound on few — > many, 



a. 



inst 



(E) < a mst E, m 



-) 
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for any v 



(8.22) 



in the regime ( |S.20 ). Furthermore, one argues that the exponent F(E/E ) for few — > many 
probability is obtained from F(E/E , v) in the limit v — > 0, i.e. with exponential accuracy one 

i^(A,_ )l. (8.23) 
aw \Eq J. 

This expectation has been confirmed by explicit perturbative calculations about the instanton 
339| . The reason for eq. ( |8.23| ) is that in the limit v — >• the overlap between the two-particle 
initial state and initial states with rij = v/aw particles is not exponentially small ||315| . 



aTZ n JE) oc exp 
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Thus, the idea of refs. ||337|, [338|1 is to evaluate F(E/E , v) semiclassically and then study the 
limit v — > to obtain the exponent for 2 — > any cross section. 
To implement this idea, one considers the quantity ||337| , |338 
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hist 



(E,rii 
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(8.24) 



The sum runs over all initial states obeying the constraints 

N\i)=m\i) 
H\i) = E\i) 
P|z) = 

where N, H and P are the particle number operator, Hamiltonian and operator of total spatial 
momentum in the Fock space of initial states. The final states in the sum (|8.24j) are arbitrary, 
and S is the ^-matrix in the instanton sector. The quantity a(E,ni) may be viewed as the 
"microcanonical" probability for instanton transitions from states with given cm. energy and 
number of incoming particles. It is this quantity that enters eqs.( |S.2lD and ( |8.22j ). 

A convenient basis in the space of initial states is provided by coherent states \a). Let us 
recall that the S'-matrix is given, in the coherent state representation, by the following functional 
integral |4(| |4| 



(b\S\a) = S(b*,a) 



DfaDcfrfDcj) exp B^fc, a) + B f ((p f , b*) + iS TttTf 



(8.25) 



where the integration is over initial (t = Tj), final (t = Tf) and intermediate values of the field(s), 
the boundary terms are 



B; 



d'k a k a-^- 2iuJkTi - I d 3 k ^(k)^(-k) + / d 3 k v^/K^-kje '~ 
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d?k u k f (k)<i) f (-k) + J d 3 k v / 2^fefe0/(-k)e 



iu) k Tf 



and the limit T, — > — oo, Tf — > +oo is assumed. All bosonic fields are denoted generically by 0. 
Summation over states is represented by the integration over the coherent state variables with 
exponential weight 



E 



Da,kDa* k exp 



d 3 k a* k dk 



The microcanonical probability can be written in the following form: 

a(E 1 n l )=J2\(b\SPp,P ni \a}\\ 

i,f 



(8.26) 
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where Pp M and P ni are projectors onto the subspace of a fixed number of incoming particles and 
a fixed total four- momentum P M = (E, 0). Summation in eq. flB.27j ) runs over all initial and final 
states. It can be shown ||342| , |337|| that the matrix elements of the projection operators are 



(a\P P \a) 



d 4 £ exp -iP^i 



d 3 k a* k ak.e 



ikt, 
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(a\P ni \a) 
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dr\ exp ( —iriiT] + / d 3 k a* k a k e 



(8.28) 



Combining eqs. ( |3.25| ), QS.26| ) and ( |3.28| ) and performing trivial integrations and changes of 
variables, one obtains the double functional integral representation for the probability, 
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The integrand here does not depend on (£ + £') due to translational invariance; the integration 
over d(£ + £ r ) produces the usual volume factor. Similarly, the integration over d(r] + 7]') produces 
an irrelevant pre-exponential factor. 

The remaining integrations are of the saddle point character in the regime ( |S.20|) : upon in- 
troducing the variables <fi = g<fi and a, b = ga, gb, all terms in the exponent become explicitly 
proportional to 1/g 2 . So, the general form of the probability, eq. ( |8.21| ), is immediate. Further- 
more, the exponent in eq. flS.21|) , {Ait / otw)F(E / Eq, u), is equal to the extremum value of the 
exponent in eq. ( |8.29| ). 

Thus, one has to extremize the functional 

$(0,0';a,a*;6,6*;T; 9) = 



ET + m9-J d 6 k ala k ^ kT+e - J d 6 kb* k b k 

+ B^, a) + B f {<l> f , b*) + B*M, a*) + B* f (<f>' f , b) + iS(4>) - iS ((/>') 
with respect to all its variables. Here we introduced the notation 
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Without loss of generality one considers only real values of T: the imaginary part of T may be 
removed by time translation. Perturbative calculations ||338|| suggest also that the saddle point 
value of 9 is real. Extremization of the functional ( 8.30 ) is conveniently performed by moving 
the contour in the complex time plane from the real axis to the contour ABCD shown in fig. O. 



At this contour one obtains the following boundary value problem ||343 
i) The field obeys the usual field equations, 



5_S_ 

50 



0. 



ii) In the future asymptotics, region D, the field is real, i.e. it is real on the line CD, 



I CD 



real. 
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iii) In the past asymptotics, t = iT/2 + t,t = real — > — oo (region A), the field is a collection 
of linear waves whose positive and negative frequency parts are related to each other 

0(k, t) = f k e^ + e- e r k e-^. (8.31) 

The values of T and 9 are related to E and rii through 

E = [ d 3 k u k f k f* k e- e (8.32) 

Hi 



J d 3 kf k f* k e- 9 . (8.33) 



Equation (|8.32|) is the natural result of the fact that the energy of the classical solution is equal to 



the energy of the scattering process, while eq. (8.33) is the analogue of the usual relation between 



the Fourier components of linear classical fields and the corresponding number of particles. The 
exponent F(E / Eq^v) is equal to the extremum value of the functional ( |8.30| ) evaluated along 
the contour ABCD (in fact, only the part ABC is relevant). Of course, one has to make sure 
that the solution to this boundary value problem indeed describes the instanton-like transition, 
i.e. that the topological numbers of initial and final vacua differ by one. In fact, an appropriate 
quantity in the case of a finite number of incoming and outgoing particles is the winding number 
of the Higgs field |3l7| . 



As expected, requirements (i)-(iii) represent a purely classical field theory problem. However, 
the fields are necessarily complex on the contour ABC: eq. ( |8.31| ) is the spatial Fourier transform 
of a complex field. Also, the field must have singularities somewhere between the real axis and 
the line AB, otherwise the conditions (ii) and (iii) would contradict each other. It is natural 
that when F(E/Eq,u) < (we assume implicitly that this is indeed the case), the classical 
problem is formulated on the contour in the time plane that contains both Minkowskian and 
Euclidean parts: we are dealing with a kind of tunnelling ( "Euclidean" ) process, but incoming 
and outgoing particles live in Minkowskian time. 

A special case of the above boundary value problem emerges at 

6 = 

when the field both at final and initial asymptotics is real. In this case the classical solution is 



real on the entire contour ABCD of fig. 13, and hence has turning points at t = and t = iT/2 
(points B and C) |g , 

T 

<9 t 0(x, t = 0) = <9t0(x, t = i—) = for all x. 

In other words, the solution is real in Euclidean space-time and periodic in Euclidean time 
with period T. This periodic instanton describes the instanton-like transition with maximum 
probability at a given energy E, and the corresponding number of incoming particles is the 
optimum number at this energy ||342|| . The maximum probability is determined by the truncated 



Euclidean action of the periodic instanton, 

aZl(E)<xexp[ET-S per (0,T)]. 
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At low energies, the periodic instanton is represented by the instanton-anti-instanton chain, 
while at E close to E sp h its fields are those of the sphaleron plus small oscillations in the 
sphaleron Minkowskian negative mode (which is the only positive mode in Euclidean time). At 
E = E sp h the maximum probability becomes unsuppressed. The number of incoming particles 
at which this occurs is equal to that produced in the decay of the sphaleron, 



W'sph 



1 

aw 



At E > E sp h the instanton-like transitions are unsuppressed at some n,, and the periodic instan- 
ton does not exist. Periodic instantons in various models have been obtained numerically for all 
energies in refs. ||344| , |345| , |346|| and analytically in ref . |p47| . 



The approach outlined above is of a fairly general nature. It enables one to study, at least in 
principle, various processes similar to instanton-like transitions induced by particle collisions. In 
fact, most extensive results obtained within this approach up to now refer to false vacuum decay 
induced by collisions of highly energetic particles in scalar theories. The reason is that solving 
classical field equations is easier in scalar theories. Also, the form of the scalar potential may be 
suitably adjusted. Since the false vacuum decay is in many respects analogous to instanton-like 
processes in gauge theories, it serves as a theoretical laboratory to probe various sets of ideas 
p0| , |35ll , |352l , PTt [353| , |515fl . 
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A particular model where the boundary value problem formulated above can be solved ana- 
lytically in a wide range of the parameters E/E and v is provided by the theory of one scalar 
field in (1 + l)-dimensional space-time described by the scalar potential [|347]| : 



2 2 2 

77T , 2 m A v z 
V{<t>) = -jP - — exp 




(8.34) 



where v 2 ^> 1 is the parameter that plays the role of the inverse coupling constant. A is another 
free parameter which is chosen to be large, A ^> 1 (but A < u). The potential has the form 
shown in fig.|l4| At large A it is quadratic almost up to <fi = v and has a steep cliff at <fi > v . The 
problem that is addressed is the decay of the metastable vacuum = induced by collisions of 
energetic particles in this vacuum. 

At low energies this process is described by the bounce ||127 |, which is the analogue of the 
instanton. There also exists an analogue of the sphaleron, which is a critical bubble [126]. Its 
energy equals the height of the barrier separating the two phases (the "true vacuum" in this 
model is = oo, but this pathological property is irrelevant for the problem of the false vacuum 
decay). The sphaleron energy in this model is 
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sph 



mv 



and the characteristic number of incoming particles is 
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-v 2 . 
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At E > E sp h and rii > n sp h the induced false vacuum decay occurs without exponential sup- 
pression. The parameters appearing in the classical boundary value problem are E/E sp h and 
n i/n sp h, and the probability of the induced vacuum decay has the form 



a(E, rii) oc exp 



SrE 
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-,z/;A 



sph 



where Sb = const • v 2 is the action of the bounce and F(E/E sp h — > 0) = — 1. The question is 
whether the exponent F becomes zero at some energy for v < 1 and, most importantly, whether 
it becomes zero at high enough energies or remains always negative in the limit v — > 0. 

The way to solve the classical boundary value problem in this model is to solve the free massive 
field equation in the region of the space-time (on the contour of fig. |13|) where <fi < v, solve the 
massless Liouville equation in the region where <p > v, and then match the solutions ||347|| . In 
this way it is possible to analyse the initial multiplicities which are not too low, 
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The result for the function F as function of energy at various v is shown in fig. [15 . 

Figure [15] shows that there exists, at rii obeying eq. (|8.35| ), some critical energy E crit {ni) at 
which F becomes equal to zero, and the exponential suppression disappears. At 1 ^> n,i/n sp h 3> 
A -1 the expression for this critical energy is fairly simple: 



p , , 4 / 7T 2 n sph 
Ecritirii) = - exp 
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Clearly, E crit rapidly grows as rii becomes small. These results show that at rii < n sph one 
can still have unsuppressed instanton-like transitions at the expense of increasing the energy. 
Furthermore, at sufficiently large nj, namely at rii obeying eq. (|8.35|) , there exist real classical 
solutions to the field equation in Minkowski space-time ||353|| that describe the false vacuum 
decay induced by rij incoming particles in purely classical (hence unsuppressed) manner at 
E > E crit {ni). 

Unfortunately, analytical solutions to the boundary value problem are lacking at ni/n sp h ^ 
A" 1 . In particular, the limit v — > cannot be studied even in the specially designed model 



However, it has been shown ||353|| that there are no classical Minkowskian solutions describing 
unsuppressed induced false vacuum decay at n,i/n sp h < n 2 /A and any energies. This means that 
the false vacuum decay induced by collisions of two particles is exponentially suppressed at all 
energies, in accord with the general unitarity argument. 

The above boundary value problem is suitable, at least in principle, for numerical calculations: 
after all, one has to solve classical field equations with specified boundary conditions. The 
numerical study of induced false vacuum decay has been undertaken in ref. ||346|| in the context 
of (3 + l)-dimensional scalar theory with the potential 



m , 2 ^ iA 
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which is similar to fig. [L4[ The main result of this study is summarized in the plot of the lines 
of constant F(e, v) in the (e, v) plane, where e = E/E sp h and v = ni/n sp h (here E sp h ~ m/X and 
n sp h ~ 1/A are again the energy of the critical bubble and the characteristic number of incoming 
particles). This plot is shown in fig.[T(3. 

The most interesting region, !/<l, was not accessible to numerical calculations for technical 
reasons. However, because of inequality ( |8.22| ) the results presented in fig. [H] show that the 
false vacuum decay induced by two-particle collisions is exponentially suppressed at least at 
E < 3E sp h. Indeed, the line of unsuppressed induced vacuum decay, F = 0, is above v ~ 0.4 
at these energies (i.e. rii should be larger than 0An sp h for the transitions to be unsuppressed). 
Further analysis [|346f reveals that the region of this suppression extends at least up to 10E sp h, 
and, most likely, to infinity. This is again in accordance with the unitarity argument. 

Another approach closely related to one discussed above is to consider real classical solutions 
to Minkowskian field equations, i.e. scattering of classical waves. To every classical solution 
that disperses into free waves at t — > ±00 one can assign the number of incoming and outgoing 
particles, both of which are naturally of order 1/aiw The probability of the scattering of these 
multiparticle states is not suppressed. At given energy one tries to minimize the number of 
incoming particles under the condition that the topological number changes by 1 (instanton-like 
transitions) or that the phase of the system changes (false vacuum decay). If the minimum 
number of incoming particles tends to zero (in units 1/ckw) as the total cm. energy approaches 
some E cr , then few — > many processes are not suppressed exponentially at E > E cr (this includes 
the more likely possibility that E cr = 00, in which case the exponential suppression of few — > 
many cross sections disappears at asymptotically high energies). In the opposite case when the 
minimum number of incoming particles needed to induce the classical transition remains larger 
than const • 1/aw, the exponential suppression of few — > many persists at all energies, but the 
actual exponent cannot be calculated by studying classical scattering. 

We have pointed out already that the two approaches (complex time and real time) nicely 
match in the (1 + l)-dimensional model (|8.34| ) and that, in this model, the number of incoming 
particles required for the instanton-like transitions to occur classically is finite in units n sp h- The 
classical real time approach is also fairly suitable for numerical analysis; considerable progress in 
this direction has been reported in ref. ||354 |, where it has also been found that the instanton-like 
transitions may proceed classicaly at somewhat below n sp h at sufficiently high energies. 

Let us stress that the idea of using many — > many transitions as an intermediate step to 
few — > many is not the only way (and, maybe, not the best way) to study the instanton-like 
processes semiclassically. A completely different — and promising — approach is the generaliza- 
tion of the Landau technique for calculation of semiclassical matrix elements [p55|| to quantum 
field theory ||307| , |356| , |357| , |358| , |359| , |360|| . It remains to be seen whether this approach will be 
able to provide new insight into electroweak B and L violation at high energies. 

To conclude, the existing calculations suggest the following overall picture of the instanton- 
like transition in high energy collisions. When the number of incoming particles is of order 
I /aw, the instanton-like processes occur at unsuppressed rates provided that the energy is 
sufficiently high. This is perfectly consistent with the calculations of the sphaleron rate at high 
temperatures: multiparticle collisions are possible in the high temperature plasma, and they are 
responsible for the high rate of electroweak B and L violation. On the other hand, if the number 



72 



of colliding particles is small (say, two), the instanton-like processes occur at exponentially low 
rates at all energies, which is in accord with unitarity argument. The actual suppression factor 
is, unfortunately, still unknown for the most interesting range of energies. While this suppression 
factor may be of limited interest for the electroweak theory (the suppression by exp(— const /aw) 
with almost any constant in the exponent will make the transitions unobservably rare), it may 
become important for the study of similar processes in QCD (for the discussion of instanton- 
induced processes in QCD see refs. |ft61|, |362|, |363|, |364|| ) . 



9 Conclusion 

The non-conservation of baryon number in the early Universe, giving rise to the baryon asym- 
metry, was proposed by Sakharov almost thirty years ago. Still, its particle physics origin is not 
established. Certainly, there were baryon number violating interactions operating at tempera- 
tures well above 100 GeV — these were anomalous electroweak reactions — but whether the 
baryon asymmetry came entirely from this source or was a combined effect of the electroweak 
processes and grand unified or/and intermediate scale interactions is unclear at the moment. In 
any case, the explanation of the existing baryon asymmetry requires at least a mild extension 
of the Minimal Standard Model. 

One expects substantial further progress, in coming years, in the understanding of the elec- 
troweak baryon number non-conservation and its role in the early Universe. On the theoretical 
side, quantitative estimates of the regions in the parameter space where the observed amount of 
baryon asymmetry is produced, are to be obtained in various extensions of the Minimal Standard 
Model. This requires further development of the kinetics of the electroweak phase transition, 
B-violating reactions and non-equilibrium description of baryogenesis. The calculations of the 
suppression factor for the anomalous reactions in high energy collisions are to be done in MSM 
and its extensions, and the relevance of similar reactions for hard processes in QCD is to be 
understood in detail. Most remarkable progress is expected, however, on the experimental side. 
Uncovering the physics at the energy scale of a few hundred GeV to a few TeV (Minimal Stan- 
dard Model? Supersymmetry? Extended Higgs sector? Technicolor-type symmetry breaking? 
...) will be crucial for establishing or ruling out the electroweak origin of the baryon asymme- 
try of the Universe. Finding out the mechanism of CP non-conservation in Nature (5-meson 
physics, electric dipole moments of neutron and electron) will become another important step. 
Possible experimental discovery of lepton number violation (Majorana neutrino masses, neutrino 
ocsillations, muon number non-conservation in /z-decays) and/or proton decay would be a strong 
indication of the early origin (T ^> 1 TeV) of the baryon asymmetry. Having understood, to a 
considerable extent, highly non-trivial aspects of the electroweak physics in the early Universe, 
one needs strong experimental input to solve one of the most challenging problems in cosmology. 
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Figure 1: Schematic plot of the static energy as function of gauge and Higgs fields. The minima 
correspond to the classical vacua. 
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Figure 3: The critical value x c = A3/ ' g\ as a function of the physical Higgs mass m# and the top 
quark mass mt op - In general, x depends on the Higgs mass, the top mass and logarithmically on 
the temperature. From ref. [|169| 
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Figure 5: The evolution of the distribution of (i?f,) with (3r (temperature). From ref. [|16S| 
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Figure 6: The probability distribution of the average Higgs length squared R 2 used for a surface 



tension determination. From ref. 169 . 
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Figure 7: The temperature evolution at the electroweak phase transition in the adiabatic case. 




Figure 8: Leading order contribution to the instanton- induced amplitude. 
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Figure 9: First corrections to the amplitude 2W — > nW. Dots represent the residues of the 
instanton field, solid lines correspond to residues of the propagator, in the instanton background, 
on the mass shell for both momenta. 



Figure 10: Higher order hard-hard correction. 




Figure 11: VF-boson propagator at low Q 2 with perturbative and non-perturbative (instanton- 
anti-instanton) contributions. 
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Figure 12: Total cross section with perturbative and instanton contributions. 
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Figure 13: The contour in complex time plane appropriate for the formulation of the boundary 
value problem for many — > many transitions. 




Figure 14: The scalar potential with unstable vacuum at = 0. 
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